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ABSTRACT 


There exists a demand and supply of goods at various 
stations at different times. There are constraints like capacity 
of the tracks, number of empty wagons available. Along with the 
cost of transportation there exists a cost associated with delays. 
Railways have to route and schedule not only goods loaded wagons 
but also empty wagons, so as to minimize costs and delays. 

A mixed integer linear programming formulation has been 
given. Literature survey is mainly related to minimum cost flow 
problems in network and vehicle routing problem. We have 
developed a heuristic for a single good wagon scheduling problem. 

The methodology consists of dividing time horizon in the 
problem into a number of periods. Now in each time period 
scheduling of goods wagons and empty wagons is done. In case of 
goods wagon scheduling we try to solve a min-cost flow problem for 
a cumulative demand and supply in that period. Then routes are 
found fitting to the solution of min-cost flow, by simulation. 
From the scheduling of the goods wagons we come to know about 
requirements or availability of empty wagons. Allocation of empty 
wagons, is done by solving transportation problem with a criteria 
of minimum shortest distance. Simulation used in the routing and 
scheduling of empty wagons, is based on choosing next available 
shortest route. Procedure is repeated to cover all the periods in 
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the time horizon. 
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CHAPTER 1 
INTRODUCTION 


1*1 INDIAN RAILWAYS 

The first railway started in India in 1855. It ran 
over a route of 54 km long, from Bombay V.T. to Thane. Since then 
Indian Railway's network has grown enormously. Today, railways 
run over a route 70,000 km. long. It is the largest railway 
network in Asia and also fourth largest in world. Operations of 
railways have grown tremendously. It has emerged as a largest 
organization in service sector in India. In the world Indian 
Railways is second largest railway under a single management. To 
carry out the operations of railways, it is required to employ Z0 
lac employees, which is more than any other organization in India. 
Even though Railways, have enormously grown, both passenger and 
goods traffic has grown at a much higher rates straining the 
resources of railways. This has imposed many constraints. 

The Indian Railways haul 40 crore tonnes of goods in a 
year. Mainly it is the freight movement, that contributes to the 
profit of Railways. In fact Indian Railway bears a loss of more 
than Rs . 1500 crores on passenger transportat ion. Hence it would 
not be wrong to say, that goods transportation 


is of primary 
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importance to the railways because of the revenue it earns and 
profits it contributes. So also the scheduling of goods trains is 
the area of interest of our thesis work. 

!• 2 CHARACTERISTICS OF THE GOODS SCHEDULING PROBLEM 

The passenger trains run according to the fixed time 
table. Even the new trains which are introduced are scheduled at 
time mostly according to the convenience of the passenger. Very 
few trains are introduced every year. As such, doing the schedul- 
ing technically, has very little scope. Contrary to this, goods 
trains do not have fixed time table, not even fixed routes, halts. 
The railways have to do entire scheduling (both day to day and 
long run) of goods trains. 

The Indian Railways is owned by Government of India. 
Hence, many a times it has to look after not only the movement of 
the goods, but also the distribution of the goods. It includes 
goods like food grain, petroleum products, fertilizers, cement, 
steel so on. Hence the nature of problem is some what tougher 
than the case in which customer would like to send goods from one 
place to another. This is because it is required to consider 
empty wagon scheduling and also time delays and accompanying cost 


i ncreases . 
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i* 3 THE OUTLINE OF THE PROBLEM 

In a typical transportation problem from operations 
research, we have warehouses which can supply certain quantity of 
goods. Also there are markets which have certain amount of demand 
to be satisfied. The cost of transporting from each warehouse to 
every market (cost matrix) is known. Similar to transportation 
problem, we have in our railway network supply nodes (stations) 
and demand nodes. We know for each arc in the network the cost of 
transporting both loaded wagon trains and also empty wagon trains. 
It should be noted that in our case cost of transportation can not 
be simply specified as cost of transportation from a supply node 
to demand node. This is because we have to consider not only the 
cost of transporting loaded goods but also that of empty wagons. 

Amongst the many constraints we have the important one 
being the capacity of the rail tracks, expressed in terms of 
number of wagons / trains that can be moved over a track in 
particular direction in unit time. This makes our problem 
partially is similar to a min-cost flow problem, having many 
source and sink nodes, instead of just one origin and one destina- 
tion node. However, the complexity of the real life problem of 
railways is of orders magnitude more than that of a min-cost flow 
problem. The railways have to deal with many types of goods which 
makes problem somewhat similar to the multi commodity min-cost 
flow problem. Further various types of goods can be transported 
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only in particular type of wagons. Each possible good type and 
wagon type combination having its own cost of transportation. The 
most important difficulty lies further when the concept of time is 
introduced. Every demand has a time deadline, i.e. the latest 
time by which demand has to be satisfied. So also supplies are 
available only after certain time and that supplies can be picked 
up any time after they become available. This concept of time is 
very important which is normally absent in of the typical trans- 
portation problem, travelling salesman problem, min-cost flow 
problem , vehicle routing problem. In our thesis work, we have 
dealt with the difficulties faced by the introduction of time 
associated with demand and supply. 

In our context, scheduling involves deciding the number 
of both goods wagons and empty wagons or trains moving from a 
origin node to a destination node, including the path along which 
they move, also deciding the starting time from the origin node, 
and finding the travelling time along the path, which involves 
specifying the idle time if any at any of the nodes on the path, 
on its way to the destination. Giving the delays in reaching the 
destination if any in reaching the destination. 

We use the entire model developed, to analyze how 
introduction of more resources like wagons can reduce time delays. 


cost etc. 
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1 • 4 ORGANIZATION OF THE THESIS 

Rest of the thesis is organized as follows. In the 
chapter 2, we give mathematical formulation of the problem. Here 
we have mixed integer linear programming formulation. This 
formulation reveals the complexity of the problem. In the chapter 
3, we have reviewed the literature on related topics. In the 
chapter 4, we begin by giving the reasons for choosing the 
particular methodologies, from among the various methodologies. 
This is followed by complete presentation of the methodology we 
have adopted. In the chapter 3. we take up a illustrative 
numerical problem. Then we present the computer results obtained, 
by running the program. The chapter 6, shows how the model 
developed can be used for deciding the critical investments like 
that of wagons. Finally, we present the conclusions derived from 
our results. Also we have suggested directions for future work in 


the last section. 



CHAPTER 2 

FORMULATION OF THE PROBLEM 


2» i INTRODUCTION 

In the previous chapter, the outline of the Railway 
Scheduling Problem was presented. In this chapter we present a 
mixed integer linear programming formulation of the problem. 

There are demands and supplies for different types of 
goods at different stations on the railway network at various 
times % We have to schedule the goods wagons and empty wagons so 
as to satisfy the demands and supplies. Decisions regarding the 
movement of wagons and the type of goods they carry are to be 
taken. Loading and unloading schedules of wagons over a period of 
time are also to be decided. 

2* 2 FORMULATION 

2» 2* 1 Notations 

Decision Variables 

WUL . . No. of Wagons of type w that get unloaded at i with 

l kwt 

goods of type k on day t. 

WL^wt No. of wagons of type w that get loaded at i with good k 
on day t . 

WLM. No. of loaded wagons of type w that start moving from i 

i jkwt 
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WLI „ . 
i kwt 

WEM . . 

1 jwt 

WEI . . 

lWt 

DEL ikt 

CL i j kw 
CE 

1 J W 


WAV 

w 

D ikt 

S ikt 

C kw 

CD., 
i k 

I 

T 

W 


to j carrying goods of type k on day t. 

No. of loaded wagons of type w that remain idle at i 
with goods k on day t. 

No. of empty wagons of type w that start moving from l 
to j on day t. 

No. of empty wagons of type w that remain idle at i on 
day t . 

Shortfall at node l, for goods of type k at time t. 

Constants CData) 

Cost of carrying from i to j the goods of type k in the 
wagons of type w. 

Cost of carrying from l to j empty wagons of type w. 

Time (days) required to move wagons from i to j . 

Maximum no. of wagons of any type that can move between 
l and j on any day. 

Total no. of wagons of type w that are available at the 
beginning and also in the entire planning horizon. 

Demand at node i, for goods of type k, at time t. 

Supply at node i, for goods of type k, at time t. 
Capacity of wagons of type w to carry goods of type k. 
Cost of delay at node i, for goods of' type k. 

No. of stations (nodes) on the railway network. 

Planning horizon in days. 

Number of different types of wagons. 
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2 . 2.2 

Objective 

Minimize 


A Linear Programming Formulation 
Function 


E £ £ E CL. . .<£ WLM. . + r r 

» i k W 1 Jkw t 1 Jkwt 7 7 


V (CE 


^ V LIPM 


E E CO <E DEL > 
i k t Kt 


Constraints 

(1) Empty wagons balance 

Sum of number of empty wagons unloaded, idle on day t-land 
arriving on day t, a node i , on a day is equal to the number of 
wagons loaded, idle, left on day t at node l. Here it is assumed 
that loading / unloading takes one day. e.g. Wagon unloaded on the 
day t-1 can not leave on the same day either as empty or as loaded 
wagon but can leave on next day or later. 


E WUL WEI ,...+£ WEM . ...... 

r* lkw(t-l) iw(t-l) . jiw(t-tjt) 

K J 

= E WL.. . + WEI. . + £ WEM. . 

*f l kwt i wt — i iwt 

k j 


for al 1 i = 1 . . I 
w = 1 . . W 
t = 1 . . T 


S imi larly , 

<2) Loaded Wagons Balance 

Sum of number of wagons loaded, idle and arriving at a node 
i, on a day is equal to the number of wagons unloaded, idle and 
left node i . 


£ WL . . + r WLI . ,, ,.+ ££ WLM . . , . . . . 

^ ikw(t-l> — lkw(t-l) " " jikw(t-tjv) 


= E WUL.. . + £ WLI .. . + £ £ . . 

7 . lkwt " lkwt . r i jkwt 


WLM 


j k 


for all i = 1 . . I 
w = 1 . .W 
t = 1 . .T 
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(5) Satisfying demand and supply 

The cumulative capacity of wagons unloaded for satisfying 

demand at any node i, should be greater than the demand existing 

for all the types of goods, at any time, 
ti ti 


£ WUL 
t = 0 


1 kwt ' 


'kw 


+ DEL . . . 
i kt 


> 


E D 
t = G 


i kt 


for demand 
for all ti = 1 

V i , k. 


T 


Similarly, the cumulative capacity of the wagons loaded at a 

node should not exceed the supplies available at that node i, for 

all types of goods, at any time. 
t± ti 


E WL 
t = 0 


i kwt 


'kw 


E S . . for supply 

t =0 for all ti = 1 ,.T 

V i , k. 


(4) Upper bound on wagons in transit 

E WEM ... + E WEM - • . + E E WLM . . . + E E WLM . 4 . 

w 1 w ’ k w 1 k w 


< U 


for all links (joining i&j ) 
for a 1 1 t 


We assume that loaded wagon which has to be unloaded, 
can not remain idle at a node (station). It has to be unloaded 
and kept idle. In other words the loaded wagon can remain idle at 
a node only if it is to be transported to some other node sooner 
or later and not if it is to be unloaded. Hence no. of unloaded 
wagons has to be less than the no. of wagons which arrive at a 
node for same type of wagon carrying same type of goods 


WUL.. . 
i kwt 


< E WLM . „ , . . .. x 

y j i kw( t-tjt ) 


for all k , w & t 
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(6) Total no. of wagons in the network is fixed 

£ £ ^*~*ikwt + ^ WEI ,wt = WAV w for w ancl t = 0 

l k i (those at starting) 

2.3 OBSERVATIONS FROM THE FORMULATION 

In the last section we have given a mixed integer linear 
programming formulation for the problem of goods and wagon 
schedul i ng . 

It may be seen that we have included cost for not 

meeting the demand in time. The costs CD could be appropriately 

i K 

chosen to reduce the delays. However, if we try to solve the 
problem according to above formulation, then we find that even for 
a small network, the number of variables and number of constraints 
is very large. For example, if we want to solve a problem which 
involves a railway network of say 10 nodes and 20 arcs, only 4 
types of goods to be transported, in 5 types of wagons and our 
time horizon is restricted to 100 discrete units of time on which 
wagon trains can be scheduled. WLM f alone requires 16,000 

variables to be defined, apart from the other decision variables 
to be defined and constant data to be given. The total number of 
constraints is around 25,000. This becomes extremely difficulty 
to solve even if we approximate all the variables as continuous 
variables rather than integer variables. This is well beyond the 
limits of the widely available linear programming packages. 
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Solomon C 29] has shown that a problem of finding a 
feasible solution to a problem of a similar nature (but with 
lesser complexity) is a NP complete problem. This issue is 
discussed in the next chapter on literature survey. We limit 
ourselves to developing a heuristic method (which is computation- 
ally viable) which will generate reasonably good schedules. Hence 
the main aim of our work, is to develop algorithm and write a 
computer program for the same. The chapter 4 is devoted to 
explaining the methodology developed. 



CHAPTER 3 

LITERATURE SURVEY 


•3* 1 INTRODUCTION 

We begin the literature survey by studing the mincost 
flow problem, where the concept of time does not come in. Later, 
we study the work of various shades of Vehicle Routing Problem. 


3.2 MINIMUM COST FLOW PROBLEM CMCFP) 

The initial description of the transportation problem is 
due to Hitchcock <1941). Since then, there have been a number of 
papers concerning the computational aspects of MCFP. Surveys of 
earlier work were done by Charnes and Cooper ( 1 951 > , Ford and 
Fulkerson [14], Dantzig C9], Busacker and Saaty C4], Hu CIS], and 
Golden and Magnanti C16]. 

One of the most important and active area of interest 
today is the basic simplex solution techniques which form the main 
tool of the work here. The approach is to use a basis tree to 
implement a procedure analogous to the primal simplex algorithm 
for the general linear programming problem. Studies have indica- 
ted that this is perhaps the most efficient way to solve pure 
minimum cost flow problems. The early precursors to this work 
were Dantzig <1954) and Charnes and Cooper (1951) who applied the 
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simplex technique to the transportation problem. In the seventies 
several authors have utilized the special representation of the 
basis trees to increase markedly the computational efficiency of 
the primal simplex techniques. Langley et al.[203 uses the triple 
label method in a primal simplex algorithm for the capacitated 
transportation problem. Glover, Karney and Klingman [153 apply 
the triple label technique to the pure network problem. They have 
also shown the superiority of the primal simplex procedure for 
network flow programming problems which was attested by others. 

Another approach, which has also been incorporated in 
this work, as an alternate solution strategy is the incremental 
flow procedure. An incremental approach is discussed in Ford 
and Fulkerson [143. Assad A. A [13 gives comprehensive survey of 
the literature dealing with multicommodity flow problems. A 
decomposition and partitioning approach is used to solve the 
min-cost flow problem using arc-chain formulation of the problem. 
The arc-chain formulation of the linear minimum cost flow problem 
suffers from the defect that the chains in the network should be 
enumerated in advance. These may be enormous in number and thus 
preclude solving a large network problem. An early paper of Ford 
and Fulkerson [153 suggested that the chains be generated as 
neeeded according to a column generation techniques. Tomlin C523 
adapts the procedure in the context of the Dantzig - Wolfe 
decompos i t l on . 
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3 ‘ 3 VEHICLE ROUTING PROBLEM 

The vehicle routing problem (VRP) involves the design of 
a set of minimum cost routes, originating and terminating at a 
central depot, for a fleet of vehicles which services a set of 
customers with known demands. Each customer is serviced exactly 
once and furthermore, all the customers must be assigned to 
vehicles such that the vehicle such that the vehicle capacities 
are not exceeded. The VRP area has been intensively studied in 
the literature, seeMagnant i C251. Chnstofides [5] has given 
optimisation algorithms and also heuristics for the basic VRP. WE 
refer to Solomon and Desrosiers C 3 1 3 have described various shades 
of VRP with Time Window Constraints. 

3*3.1 Vehicle Routing Problem with Time Windows 

Solomon C 3 1 3 has given mathematical formulation in which 
VRPTW and many other problem problem variations are derived from 
pickup and delivery problem (see section ?.>.£). This formulation 
appears in the Appendix A. In the VRP with time windows (VRPTW), 
the above issues have to be dealt with under the added complexity 
of allowable delivery times, or time windows, stemming from the 
fact that some customers impose service deadlines and earliest 
service time constraints. In these problems, the spatial aspect 
of routing is blended with the temporal aspect of scheduling, 
which must be performed to ensure the satisfaction of the time 
window constraints. The service of a customer, can begin 



within the latest time when the customer will permit the start of 
service. The times at which services begin are decision vari- 
abl es . 

It is worth distinguishing now between hard and soft 
time windows. In the hard time window case, if a vehicle arrives 
at a customer too early, it will wait. In addition, due dates 
cannot be violated. In contrast, in the soft time window case, 
the time window constraint can be violated at a cost, in a manner 
similar to dualizing constraints in Lagrangian relaxation context. 

The vehicle fleet size and mix, may assumed homoge- 
neous. However, many of the methods proposed can be extended to 
consider a heterogeneous fleet size. Furthermore, the number of 
vehicles used can be free, i.e., the fleet size is determined 
simultaneously with the best set of routes and schedules rather 
than being fixed a priori. 

Additional complexities encountered in the VRPTW are 
length of route constraints arising from depot time window. It 
may be noted that the vehicle departure times from the depot are 
decision variables. Furthermore, precedence relationships among 
certain customer are encountered in some problems, such as those 
involving both pickup and deliveries. 

Most of the effort has been directed at the operational 
problem of determining the best set of routes and schedules. In 
the presence of time window, the total routing and scheduling 
costs include not only the total travel distance and time costs 
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considered for routing problems, but also the cost of waiting 
time which is incurred when a vehicle arrives too early at a 
customer location or when the vehicle is loaded or unloaded. 
Several authors have also have analyzed the strategic question of 
the optimal fleet size. Keaton C 2 1 3 dealt with designing optimal 
railroad operating plans. Various related decisions are modeled 
as mixed integer linear program. Solutions are obtained using 
Lagrangian Relaxation and Heuristic approaches. 

The VRPTW has emerged as an important area for progress 
in handling realistic complications and generalizations of the 
basic routing model (Schrage [273>. In the last few years we have 
witnessed the development of a fast growing body of research 
focused on the vehicle routing and scheduling problem structures 
with time windows. Given the intrinsic difficulty of this problem 
class, the later work on the generic VRPTW has typically focused 
on heuristic methods. The findings of Solomon [293 indicate that 
this problem class is fundamentally more difficult than the VRP. 
The VRPTW is NP-hard (by reduction from the VRP) and heuristics 
have so far offered the most promise for solving realistic size 
problems . 

Solomon [29,?03 has designed and analyzed a variety of 
route construction heuristics for the VRPTW. The results of the 
extensive computational study reported in Solomon [293 indicate 
that a sequential time-space insertion algorithm proved to be 
very successful in a number of important VRPTW environments. Very 
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good initial solutions to a large number of 100 customer test 
problems were obtained. The excellent performance of such 
heuristics can be explained by realizing that while routing 
problems seem to be driven by the assignment of customers to 
vehicles as indicated by the success of the Fisher and Jaikumar 
[12], generalized assignment heuristic- the sequencing aspect of 
the problem seems to drive routing problems dominated by time 
windows. It is this aspect of the problem that insertion approa- 
ches capture so well. A very large computational requirement was 
also reported by Cook and Russell [8], in early work on the VRPTW 
for a k-optimal improvement heuristic, which was effective in 
solving an actual problem with a few time constrained customers. 
Analytical results concerning the behaviour of VRPTW approximation 
methods are derived through worst-case analysis of heuristics in 
solomon [29D. For a variety of heuristics, including improvement 
methods, it is shown that their worst-case behaviour on n-customer 
problems, is at least order of n for minimizing the number of 
vehicles used, the total distance travelled and the total schedule 
time. While heuristics have been found to be very effective and 
efficient in solving a wide range of practical size VRPTW, few 
optimal approaches have been reported, see Kolen et al. [22]. The 
former authors extend the shortest q-path relaxation algorithms of 
Christofides et al. [6] to the problem with time window. The 
largest problem solved to optimality involved 4 vehicles servicing 
14 customers with tight time windows. 
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Fer land and Fortin [113 have given a heuristic approach 
for the vehicle scheduling problem with sliding time windows The 
advantage has been taken of the fact that the starting time of 
each task must fall within a given time interval rather than being 
fixed. The approach consists of identifying pairs of tasks 
offering good opportunity costs for reducing the overall cost and 
finding the ways of modifying starting times to allow them to be 
1 inked. 

3.3.2 Pickup and Delivery Problems CPDPTVD 

PDPTW is a generalization of the VRPTW which is concern- 
ed with the construction of optimal routes to satisfy transporta- 
tion requests. each requiring both pick and delivery under 
precedence, capacity and time window constraints. (Mote that the 
VRPTW is the particular case of PDPTW where the destinations are 
all the common depot. 

Most of the literature on pick-up and delivery problems 
has appeared for the dial-a-ride problem (DARP) . It can be seen 
that the single vehicle DARPTW is a constrained version of the 
classical travelling salesman problem. Psaraftis [24] presents 
two dynamic programming algorithms using respectively backward and 
forward recursions. The objective function minimizes the total 
customer inconvenience and "maximum position shifts” define time 
window constraints. These algorithms require 0(n 2 ? n ) time (where 
n is the number of customers), a fact which limits the tractable 


problem size to no more than 8-10 customers. 
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Sexton and Bodin [28] give mixed integer nonlinear 
programming formulation for the single vehicle dial-a-ride problem 
with given maximum delivery times or specified minimum pickup 
times to minimize customer inconvenience. Benders' decomposition 
is used to partition the problem into a routing component (a 
linear program). For the scheduling problem, an optimal algorithm 
which exploits the network flow structure of its dual is develop- 
ed. The overall procedure is heuristic due to the difficulty of 
solving the routing problem exactly. Its computational tractabi- 
lity stems from the non-iterative nature of the scheduling. 
Successful computational experience on moderately sized real data 
is reported. Problems with 7 t - o 20 customers were solved in an 
average of 18 seconds of UNIVAC 1100/81A CPU time. 

Jaw et al . [19] also use a parallel insertion algorithm 
which appears to be very effective and efficient in minimizing a 
weighted combination of customer dis-utility and system cost. The 
authors report computational experience with simulated and real 
data on the VAX 11/750. The former problems. involving 250 
customers and 10-14 vehicles, took around 20 seconds each. 

The Bodin and Sexton [51 procedure is a traditional 
"cluster first, route second" approach. For a fixed fleet size, 
it partitions the set of requests into vehicle clusters and solves 
the resulting single vehicle dial-a-ride problems using the 
heuristic based on Benders decomposition (Sexton and Bodin C281). 
Requests are then moved one at a time from one vehicle to another 
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while attempt i ng to reduce total user i nconven 1 ence . 

3.3.3 The Multi-Period Vehicle Routing Problem CMPVRP3 

Here, the time windows are full days and a service 
activity must occur on a specified number of days of the planning 
VRPTW where each customer has to be visited a given number of 
times within its multiple time windows. Most of the existing 
literature is based on application that require the solution of 
large scale problems. 

Beltrami and Bodin [23 address this problem in the 
context of routing hoist compactor trucks. The problem, which 
involves customers requiring service either 5 or 6 times per week, 
is solved heur i st l cal 1 y . One approach involves developing the 
routes first and then assigning to days of the week. In the 
second approach, routing is performed after the customers have 
been randomly pre assigned to days of the week. 

Rusell and Igo [263 extened several routing heuristics 
to assign accounts to days of the week such that the weekly travel 
distance is minimized. The refuse collection problem considered 
involves the routing 4 trucks through 490 location requiring 
service from 1 to 6 times per week. The original problem is 
transformed into a pure routing problem by creating multiple 
copies of each account that requires service several times a week, 
resulting in a 776 customer problem. The best approach found was 
to assign the locations to days of the week using a simple 
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heuristic, solve the resulting problem, where the spacing condi- 
tions are enforced and improve the solution through the use of a k 
1 nterchange heur i s t i c (M-Tour , Russell [2?3>. The use of M-Tour 
proved once again to be computationally prohibitive. Christofides 
and Beasley [73 develop new heuristics for the MPVRP based on 
median problem and traveling salesman problem relaxations. The 
algorithms are characterized by the use of an interchange proce- 
dure which improves on initial choices of delivery days combina- 
tions for customers. 

Haghani A. E. [173 presents the formulation and solution 
of a combined train routing and makeup and empty car distribution 
model. The model results in a large scale mixed integer programm- 
ing problem with nonlinear objective function. A heuristic 
decomposition technique was developed to solve the model. 



CHAPTER 4 
METHODOLOGY 


1 INTRODUCTION 

The problem as formulated in chapter 2 is fairly complex 
and hence we begin solving it by relaxing a class of constraints 
at a time and solve the relaxed problem by known methods which are 
efficient. We present the solution method for a single commodity 
problem, but it can be extended to the multicommodity case. 

The detailed steps in the algorithm for solving Railway Wagon 
Scheduling Problem are as given below. 

A. 2 INPUT 

(i) Network is specified in the form of arcs. 

Data about the Arc 

(a> Nodes : nodel , node2 which are joined by the arc 
<b) Distance. 

Cost for carrying loaded wagons cost for carrying empty 
wagons is assumed proportional to distance and both 
costs are computed. 

Cost is expressed as cost/rake for travelling arc. 


Assumed to be same in both directions. 
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(ii> Demand / Supply of goods (expressed in wagons (trains)) 

Input : Node, demand / supply, time . 

(iii) Length of time horizon 

length of one time period within a time horizon. 

We partition the problem into a few distinct time periods 
say equal to *n*. This can be done by a suitable heuristic 
or by simple inspection of data. We have done it by 
inspection only. 

(tv) Availability of empty wagons at various nodes in time period 
one (or in any other period) 

Input : Node, availability of wagons, time. 

4.3 GENERATE THE INFORMATION FOR GIVEN NETWORK 

For each arc, there is similar arc in reverse direction. 
Hence, ours is und i rect i onal network for which we have following 
relations holding true. 

Number of arcs leaving a node i 
= Number of arcs arriving at a node 1 . 

= Degree of a node i / 2 

= difference between pointCi] and pointCi+1]. 

(1) Arc numbers of all the arcs leaving a particular node i are 
stored in the form of a forward star network representation. 

Arcs' feaving node i will be stored in position from pointCi] to 
pointCi+13-1 in the array ftraceC]. 
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< 2 ) Arc numbers of ail the arcs arriving a particular node i 
stored in the form of a reverse star network representation. 
Arcs entering node 1 will be stored in position from pointCi] 
pointCi+13-1 m the array rtraceC]- 


4.4 SCHEDULING OF GOODS WAGONS 

Out of these *n' distinct period problem, we pick up 
first period problem and totally drop the time constraints 
supply and demand. The min-cost flow problem PI is formulated 
below. 

Minimize £ £ C. . .X. . 

<1> Flow conservation at all nodes i 

Outflow - Inflow = Supply (Demand) 

' B ‘ v ‘ 

<2> Capacity restriction on flows over the arcs 

X . . < K . V arcs 

i J * 1 


(?) Non negativity of flow constraint 

X. . >0 V ( i , j ) 

» J 

Where 

X. . Flow from node i to node j 
» J 

C Cost of transportation unit quantity from node 
i ] 

B. Supply (or Demand) at node i 

K.j Capacity of arc from node i to node j 


to node j 


are 

to 


the 

on 

as 
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The solution to the above problem gives the flows over 
the arcs in a period. Various fractions of the flow over the same 
arc can be a part of the different routes from various supply 
nodes to different demand nodes. Solve the min-cost flow problem 
as given below. 

4-. 4.1 Solve the min— cost flow for the current period 

Find the cumulative demand / supply of goods within the 
current time period. This is done because demand / supply may 
occur at different times in same time period. The demands must be 
satisfied within the same time period only and also at time 
mentioned. However supplies not exhausted during past time 
periods are considered available in the current time period also. 

Cumulative demand = 2 Demand for current period. 

Cumulative supply = £ Supply for current period + 

residual supply of previous periods 
With this information about cumulative demand / supply at each 
node and the data about the network, solve the min-cost flow 
problem for the current time period. The capacities of the arcs 
used in solving min-cost flow problem are specified for entire 
period. The min-cost flow is formulated as a linear program 
and solved by a revised simplex method. 

The output of the min-cost flow problem is in the form 
of flows over various arcs. The flow may exist in both the 
directions on a arc (arcs have been separately identified for each 


direction). Since the individual arc flows have already been 
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decided the cost of transportation is same irrespective of the 
route chosen and the decisions about which supply node supplies 
towhich demand node, since it fits in solution given by the 
min-cost flow. 

4*4.2 Solving t-he Transportation Problem Tor goods wagons 

Take a demand node for all those nodes having a demand 
and each time it has a demand. Similarly have a supply node for 
all those nodes having a supply and each time it has a supply. 
Thus more than one copies of same node will exist if and only if 
demand / supply exists more than once till end of current time 
peri od . 

Now we wish to know which node is to supply which demand 
node and also the likely quantity supplied. While doing so it is 
desirable to minimize the delays. We solve the following problem 
PJ which is the transformed transportation problem. 

Formulation of transportation problem is given here. 

Problem P? 

Maximize EES X 

** 44 i J il 

Subject to 

<1> Demand restriction for node l 

£ . > A. V i = 1 , 2 .... m 

(2) Supply restriction from node j 

E Xj j £ B V j = 1 . 2 n 

(5) Non negativity of flow constraint 

x u > 0 V <1,J> 



Compute a shortest distance matrix 
between each pair of supply node 'i 
it into a shortest time matrix ( t 


( d ) giving shortest 
i J 

• and demand node ' j ’ * 
> . 


d i stance 
Convert 


i J 

Compute a available time matrix (A. .) giving difference of time 

at which supply becomes available and the time at which demand 

needs to be satisfied. 

Compute a slackness matrix. (S ) 

i ) 

Slackness = Available time - time for travelling shortest distance 


Use the criteria of maximimum of slackness to decide allocation 
between demand and supply nodes. This criteria is used so as to 
minimize time 1 nf eas 1 bi 1 l t 1 es in current and future periods. A 
simple heuristic approach is used to solve this transportation 
probl em. 

Solution procedure to solve objective function by maximum 
slackness criteria is given below. 

Heuristic One 

Step 1 Find out slackness matrix (also known as transportation 
table) 

Step £ Select a variable with maximum slackness as basic 

variable. It is assigned a value equal to minimum of 
supply or demand. Update the residual demand and supply 
by subtracting the assigned value. 

If there is tie in the selection of the variable with a 


Step 3 
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maximum slackness, then find the penalties for not 
making that variable as a basic variable. Break the tie 
by choosing the variable which has maximum penal ty . 

Step 4 Repeat the procedure till all demand is satisfied. 

We try to stick to the above solution of transportation 
problem P3 and find the routes from supply node to demand node 
within framework of solution given by min-cost flow. 


If we get a infeasible solution to the min-cost flow 
then we need to solve a different transportation problem which 
will reasonabally take care of the costs of transportation. We 

solve a transportation problem P4 , which minimizes the sum of the 
shortest distances. 

Formulation of transportation problem which minimizes the 
sum of shortest distances is given here. 

Problem P4 

Minimize E E D X . 

La 1* , j , j 

Subject to 

<1> Demand restriction for node i 


> A. V i = 1 ,2, . . ,m 

(2) Supply restriction from node j 

£ X . < B . V j = 1 , 2 .... n 

(5) Non negativity of flow constraint 

X. >0 V ( i , j > 

i J 

Where a shortest distance matrix (D^) giving shortest 


d i stance 
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between each pair of supply node ' i* and demand node ’ j'. 

Use the criteria of minimum of shortest distance to decide 
allocation between demand and supply nodes. A simple heuristic 
approach is used to solve this transportation problem. 

Solution procedure to solve objective function by minimum 
shortest distance criteria. 

Heuristic Two 

Step 1 Find out shortest distance matrix ( transportation 
table ) 

Step 2 Select a variable with minimum shortest distance as 
basic variable. It is assigned a value equal to minimum 
of supply or demand. Update the residual demand and 
supply by subtracting the assigned value. 

Step 5 If there is tie in the selection of the variable with a 
minimum shortest distance, then find the penalties for 
not making that variable as a basic variable. Break the 
tie by choosing the variable which has maximum penalty. 
Step 4 Repeat the procedure till all demand is satisfied. 

4.4*3 Simulating to find the routes for goods wagon scheduling 

We start from a supply node and travel along any of the 
arcs leaving the node if that arc has a flow indicated by min-cost 
solution. We recursively continue search till we reach a demand 
node which has some positive allocation for that supply node and 
demand pair. When a path is found successfully, we send a flow 
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along the entire path. Quantity of flow that can be sent is 
minimum of 1. mini mum arc capacity a 1 ong the path 2. mini mum of 
supply and demand. In case we fail to reach such a demand node, 
then we backtrack to the previous node along the route so as to 
try with some other arc leaving the previous node. This back- 
tracking to the previous node along the route may have to be done 
recursively till we do not have any node left (other than supply 
node itself) in the route we are trying. In that case we give up 
the current supply node, which probably may have surplus supplies 
which may remain even after satisfying all demand. Alternatively 
its supplies are available too late and can not satisfy demand 
without delays hence considered later. 

The above procedure is repeated so as to find routes 
from other supply node to demand nodes. If the demand remains 
unsatisfied even after trying with all the supply nodes once , then 
we reject the tranportat i on solution which we developed earlier. 
It may be noted here that transportation solution developed 
earlier with criteria of maximizing the slackness was just to 
minimize any time i nf eas l bi 1 i ty that may occur during goods wagon 
scheduling or consequent empty wagon scheduling. Hence even if we 


do not stick 

to 

the allocation, it 

l s 

always 

possible to 

find 

routes from 

the 

supply node to the 

demand 

nodes within 

the 

framework of 

the 

min-cost solution 

1 . 

e. without exceeding 

the 

flows indicated 

by min-cost solution 

wh ich, in 

turn do not exceed 


the capacity of the arc specified for the period. 



If our problem is too tight problem then the infeasibi- 
lity is found while solving the min-cost flow itself. In such 

cases we use minimizing the sum of the shortest distance between 
the supply and demand node pair as the criteria for deciding the 
allocations. Now the shortest available routes are found by 
simulation. Then, this approach becomes very much similar to 
empty wagon scheduling described in details later on. However in 
case of goods wagon, their scheduling is always done so as to 
reach the destination just in time (whenever possible). On the 
other hand in case of empty wagon, scheduling is done so as to 
start from the source node as soon as they become available. 

When sending the flow along a route found 

a. Check time feasibility if any, and find the amount of delays 

when it occurs. 

b. Keep the account of residual capacity of each arc on each day. 

Goods Flows are scheduled , so as to reach demand node 
just in time of need, rather than sending the flow at a time when 
supplies are available. 

This reduces the tightness in scheduling and time 
infeasibility that might be encountered during empty wagon 
scheduling. The exception to this scheduling rule, is the 
situation when supplies are not available within required time. 
Then the goods trains are scheduled to leave the supply node only 
when supplies become available, even though they will reach demand 
node with delay. 
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EMPTY WAGON SCHEDULING 


The demand and supply of goods (goods trains) acts as a 
source of supply and demand for empty wagons respectively. 
However the time of demand / supply need not be same as that in 
data. Delays might have occurred while satisfying demand of goods 
(which leads to time 1 nf eas 1 b i 1 1 ty ) . It is possible to have a 
situation in which demand gets partially satisfied in time and 
partially with delays. Also supply of goods may be picked up at a 
later time, rather than the time at which it becomes available. 
So also a smaller portion of the total supplies available at a 
node may be picked up at different times. 

Thus to get precise information for requirement and 
availability of empty wagons, it is necessary to keep account of 
(1> Number of empty wagons and (2) time at which they are required 
or are available as soon as loading or unloading of goods takes 
place during goods wagon scheduling. 

4.5.1 Solving the tranportation problem for Empty wagons 

Similar to goods wagon scheduling we shall form demand 
nodes and supply nodes for empty wagons each time there is 
requirement for or availability of empty wagons. Then we form a 
transportation problem to decide which supply node will supply 
which demand node. However unlike previous case the basis for 
allocations is simply shortest distance matrix which gives the 
shortest distance between every pair of supply and demand node. 





The reason why we have used criteria of maximizing the 
slackness while goods wagon scheduling and that of minimizing the 
sum of shortest distance while finding allocations during empty 
wagon scheduling has to be noted. During goods wagon scheduling 
min-cost flow was solved first and the purpose of solving the 
transportation problem with maximum slackness was to minimize the 
time 1 nf eas i b 1 1 1 t i es that may occur during goods wagon and empty 
scheduling. Overall cost factor was already taken care of when we 
solved min-cost flow problem. The reason why shortest distance 
matrix is used as a basis for allocations during empty wagon 
scheduling is that we would like to minimize the costs of trans- 
portation for empty wagon scheduling also. We have also found 
that maximizing the slackness is a poor criteria during empty 
scheduling. Also it is not needed because nothing further depends 
on empty wagon scheduling, unlike the case of goods wagon schedul- 
ing from which empty wagon scheduling is derived. 

4.5*2 Deciding -the routes for empty wagons by simulation 

Procedure for deciding the routes from supply node to 
demand node for empty wagons is somewhat simpler in case of empty 
wagon scheduling than goods wagon scheduling. Here we strictly 
stick to the allocations. Moreover unlike the min-cost flow of 
the goods problem, the capacities of arc are not expressed for 
whole period but for a each unit time within a period, which 
is a fraction (to be specified) of the capacity for the period. 



We take a supply node. We try to send empty wagons 
demand node along a shortest path. In case, we find that there 
no capacity available on a particular arc in the path on 
particular day then 

(1) We try, keeping the empty wagons idle just before the arc 

which there is insufficient capacity available on that particul 

day till, sufficient capacity becomes available. 
i 2) In case we find that by keeping the empty wagons idle, be 
that particular arc, it is not possible to reach the destinat 
in time, then we try to find a next shortest path excluding 
arc(s) whose capacity has saturated. In case we encounter simil 
problem in the next shortest path also then we repeat the proc 
dure. If the distance involved in travelling by next shorte 
available path is such that it will not be possible to rea 

destination then we accept that as time infeasibility, 

wagons are scheduled along a first shortest path only. We note 
delays that are bound to occur when we send the goods by 
shortest path only. 

The steps 4.4 and 4.5 are repeated for all the time periods in 
time horizon. The concise picture of the algorithm appears 


figure 4.1 




PIG. 4.1 CONCISE PICTURE OP A LG OH I TOM 
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° OBSER V AT JONS 

The choice of length of the time period seems to be 
important. If the period length is small then the number of 
periods in the horizon increases. and there is corresponding 
increase in the computational time. If the period length is taken 
to be very large then the wagons remain idle in a period. This is 
because wagons are scheduled only twice within a time period. 
Once goods carrying wagons are scheduled and once empty wagons are 
scheduled, so as to meet the requirements for the goods wagon. 
Hence there is poor utilization of the wagons when period length 
is large. Also it may result in quite suboptimal and at times 
infeasible solutions. In which case we have try again with a 
smaller length of the period. 



CHAPTER 5 
RESULTS 


5.1 INTRODUCTION 

In the previous chapter, we have presented the solution 
methodology explaining the details of various steps involved in 
the algorithm coded. The program code has been attached in the 
appendix B. In this chapter we take up a numerical example. In 

the section 3.2, the data of the problem is given. Figure 3.1 
shows the railway network which is used for the problem. In the 
section 5 . 5 , the output results have been tabulated. In the last 
section, we have discussion on the results we have obtained. 

5.2 NUMERICAL PROBLEM 

5*2.1 Data about arcs of the railway network 


Table 5. 1 Arcs of the network 


Arc 

No . 

Nodel 

Node2 

D i stance 

( km) 

Capac i ty 

Problems 

1 2 

1 

Delhi <1> 

Jhans i 

(2) 

400 

16 

26 


Delhi <1> 

Kanpur 

( 3 ) 

400 

20 

23 


Kanpur <3) 

Jhans i 

( 4 > 

200 

3 

3 


Ahmad. (13) 

Jhans i 

(2) 

800 

5 

3 


Kanpur (3> 

Lucknow 

(4) 

100 

23 

33 
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6 

Jhans l (2) 

Bhopal (19) 

300 

30 

33 

7 

Ahmad . < 1 3 > 

Bhopal (19) 

330 

6 

6 

8 

Bhopal <19) 

Indore (20) 

200 

23 

23 

9 

Lucknow <4> 

Calcutta (3) 

1000 

1 3 

23 

10 

Bhopal <19) 

Bhusaval (13) 

400 

7 

8 

11 

Bhopal (19) 

Nagpur (16) 

330 

9 

7 

12 

Cal cut . ( 3 ) 

Jamshed . ( 6 ) 

230 

13 

23 

13 

Bombay (12) 

Manmad (14) 

230 

23 

19 

14 

Manmad (14) 

Bhusav. (13) 

200 

16 

16 

13 

Bhusa. (13) 

Nagpur (16) 

400 

7 

7 

16 

Nagpur (16) 

Jamshed. (6) 

230 

23 

23 

17 

Bombay (12) 

Pune (11) 

200 

6 

6 

18 

Manmad (14) 

Daund (10) 

230 

7 

7 

19 

Nagpur (16) 

Balhar. (17) 

200 

19 

16 

20 

Daund (10) 

Bangal . ( 9 ) 

900 

6 

6 

21 

Balhar. (17) 

Hydrab. (18) 

330 

9 

9 

22 

Jamshed . ( 6 ) 

Vishakh. (7) 

800 

6 

6 

23 

Banga 1.(9) 

Hydrab. (18) 

700 

6 

6 

24 

Bangal . ( 9 > 

Madras (8) 

330 

23 

23 

23 

Hydrab .(18) 

Vishakh. (7) 

600 

23 

21 

26 

Hydrab. (18) 

Madras (8) 

700 

3 

3 

27 

Madras (8) 

Vishakh. (7) 

830 

3 

3 

28 

Pune (11) 

Daund (10) 

100 

3 

3 

29 

Ahmad . (13) 

Bombay (12) 

300 

23 

23 



i, j : nodes 
r : arc number 


RAILWAY NETWORK USED IN THE PROBLEM . 
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Network shown contains 20 nodes (n) and 29 arcs (m) . 
There are also 29 more arcs in direction just opposite to the one 
shown in the network. These are identified by numbers from m+1 to 
2m. The cost of the wagon movement along a arc is taken propor- 
tional to the distance. 

The cost of loaded wagon train movement = distance (km) / 10. 

The cost of empty wagon train movement = distance (km) /15. 

Also to have the consistancy in the un i tsf o 1 1 ow l ng conversions is 
used . 

1 unit time = 4 hrs = 200 km. 

1 day = 6 units of time. 


5.2.2 Data about the availability of supply 

Table 5.2 Supply at various nodes 


Supply node 

Supply, Time(day) 

2 

Jhans i 

5,2; 10,5; 15,5; 4,6; 6,7; 7,8; 10,10; 5,12; 18,15 

18 

Hydrabad 

10,2; 5.4; 18,5; 10,6; 5,8; 18,9; 4,10; 6,11; 7,12 

4 

Lucknow 

2,1; 5,4; 8.5; 10,6; 12,9; 2.10; 5.12; 8,15 

9 

Bangalore 

4,1; 6,2; 7.5; 8,6; 5.8; 10,10; 10,11; 10,12; 15,15 

6 

Jamshed. 

8,2; 5,4; 10,6; 10,7; 10,8; 15,9; 8,10; 5.12 

1 

Delhi 

10,1; 12,4; 2,6; 5.8; 8.9; 10,10; 12,15 

12 

Bombay 

10,1; 10,2; 10,4; 15,5; 5,6; 10,7; 15,9; 5,10; 



10,11 ; 15,15 
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5*2*3 Data about, t-h© demand requirements 

Table 5.3 Demand at various nodes 


Demand node 

Demand, time (day) 

19 

Bhopal 

10,?; 1 ? , 6 ; 10,9; ?,10; 4.11; 8,15; 7.14 

15 

Ahmadabad 

2.6; ? , 7 ; 10,9; ?,11 ; 8,15; 7,1? 

7 

V i jaywada 

20,6; 10,8; ? , 9 ; 1?,12; 1?,15; 10.14 

8 

Madras 

1 ? , 7 ; 10,8; 20,9; 10,11; 10.15; 1?,14 

14 

Manmad 

? , 4 ; 1 ? , 8 ; 2,9; 5,10; 2,15; ?,14 

16 

Nagpur 

8,?; 7,6; 8,9; 7,10; ?,15; 10,16 

? 

Calcutta 

10,4; ? , 6 ; ?,9; 10,12; ?,15; 15,16 

20 

Indore 

5,4; 10,7; 15,9; 10,10; 10,15; 5,14; 4,16 

5 

Kanpur 

10,4; ?,?; 4,6; 10,9; 15,10; 10,15; 5.1? 


The empty wagons made available in the entire network in 
the second time period is varied for both problems. Solution is 
obtained for four cases of each problem. 

1. 17? 2. 200 5. 22? 4. 2?0 

5. 3 OUTPUT 

The output solution is volumnious. We get following 
different outputs. 

1. Distribution of number of trains delayed vs. amount of delay 
Percentage frequency of delays vs. amount of delays. 
Percentage of delay vs. amount of delays 


Average delay. 
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2. Routes followed by the trains and schedule of the trains 
seperately for both goods movement and empty wagon movement. 
It is expressed as nodes on the path and the starting time 
from each node on the path which is stated taking intc 
consideration idle time if any at a node. 

3. Load on each arc, each day. 

The maximum load on the arc, and the time at which it occurs 

4. The cost of scheduling for both goods trains and empty trains 
in each time period. And total cost of the schedule 
generated . 

5. 4 FINAL RESULTS 


Table 5.4 
Problem 1 


Number of 

Total Cost 

Average 

Empty Wagons 

of Schedu 1 i ng 

Delays 

173 

33,886 

3.13 

200 

33.613 

3.23 

223 

33.096 

4.90 

230 

33.937 

4.35 
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Table 5* 5 
Probl em2 


Number of 

Total Cost 

Average 

Empty Wagons 

of Schedu 1 1 ng 

De lays 

175 

38,786 

3.49 

200 

38,416 

3.37 

225 

38,020 

3.30 

250 

36,706 

5.01 


This results are plotted in the figure 5.2 and figure 5.3 









Total 
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s * 5 DISCUSSION OF THE RESULTS OBTAINED 

Fig. 5.2. gives a graph of total cost of transportation 
vs. number of empty wagons available. One can very easily notice 
that the total cost of the transportation decreases with the 
increase in the number of empty wagons available. It is simple to 
interpret it. When the number of empty wagons available is 
inadequate, then these need fetched from far off places before 
these can be used for goods transportation, which increases the 
total cost of the transportation. 

Figure 5 . 3 ls a graph of Delays in the movement of the 
wagons. In this graph also it can observed that there is decrease 
in delays with increase in the availability of the wagons. Very 
similar explanation seems to be valid for this behaviour also, 
i.e. more the number of empty wagons available, lesser will be the 
need to bring them from over a long distance. However a tradeoff 
is possible for the decision maker between higher cost of more 
trains and the lesser delays and transportation costs. 

If we compare the results obtained for two problems then 
we find that the problem 2 in which we incur more total cost of 
transportation has correspondingly lesser amount of delays taking 
place and vice versa. However in problem 1 we find that as the 
number of empty wagons available is decreased' from 200 to 175 
there is a slight decrease in the delays (even though the cost of 
transportation has increased). This slight discrepancy probably 



due to the fact that our 


"single pass" 


n nature. 


heuristic is 



CHAPTER 6 
CONCLUSION 
AND 

DIRECTIONS FOR FUTURE WORK 


6. 1 CONCLUSIONS 

The computational experience with the algorithm that has 
been developed seems to be reasonably well. The moderately sized 
problem, has been solved on HP 9000, in approximately 30 secs of 
the CPU time. 

We feel that the graphs that are obtained, can be very 
useful for Railways for taking decisions regarding the investments 
to be made in the critical assets like wagons. In reality usually 
there is some costs associated with the delays e.g. Overtime paid 
to the employees, more fuel and other operating costs for engines. 
If the cost of delays can be quantified, then it should be 
possible to find the sum of cost of transportation and cost of 
delays. It will help Railways for getting a idea of approximate 
investments desired so as to have the minimum net cost of trans- 


portat ion. 



48 


6* 2 DIRECTIONS FOR FUTURE WORK 

After working on the problem of Railway Wagon Scheduling 
we have realized various dimensions of this area. There is 
scope for continuing work on this problem area. 

It will be useful to work in this area to develop a 
methodology for dealing with multicommodity problem. It may be 
enough to deal with few goods. This is because 80% of the total 
freight movement is due to few types of goods like steel, food 
grain, fertilizers, cement. Attempt could be made by modifying 
methodology so as to incorporate multicommodity min-cost flow and 
multidimensional assignment problem. 

* 

Given the difficult nature of the problem it is felt 
that heuristics need to be developed. Also heuristic must be of 
the "multipass" type. So that each time we are able to improve on 
the solution of the previous iterations. 

The scope of the problem can be widened so as to a 
decision support system for other important investment areas of 
Railways e.g. Rails of the Railway lines. Certain rail routes are 
heavily loaded, where as others are scarcely utilized. The study 
of this aspect will be useful in taking decisions regarding 
increasing the capacity of the routes. 

It will be desirable to take a real life data. Also, an 
attempt needs to be made to solve at least one or two problems 
optimally. It will give us the yardstick to compare the 


perform- 



4 9 


ance of the heuristic approaches developed for such problem. 
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APPENDIX A 


Here we give a formulation for pick and delivery problem 
with time windows from Solomon C51D. Then this model is used to 
develop a taxonomy of routing and scheduling problems with time 
windows . 


Let there be n customers indexed by 1 . Associate to the 
pickup location of customer i a node i and to his delivery 
location a node n+i. Also associate to the depot node 0 and 2 n + 
1. This creates a clear distinction between customers, their 
associated locations and the nodes of the network. Note, however, 
that different nodes may correspond to the same physical location. 
Therefore. N = Cl,2,..,n, n+1 , n+2 . . . , 2n + l > is the node set for our 
network, and P = P + U P , P + = Cl,2,..,n>, P = Cn + 1 , nt2,..,2n} is 

the set of nodes other than the depot nodes. We will index both N 
and P by u and w. Customer i demands that d.^ units be shipped from 
node i to node n+i . Next, let Ca^.b^l denote the pickup tie 
window for customer i and let ta ,b D denote his delivery time 

window. Let V = £1.2 I V 1 > be the set of vehicles to be routed 

and scheduled; index this set by v. Let D be the capacity of each 
vehicle. Let also Ca ,b 3 denote the vehicles’ dparture time 

O O 

window from the depot and C a 2n+i> their time window for 

arrival back at the depot. 
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For each distinct u,w in N, let t and c represent the 

U V uv 

travel time and the travel cost form u to w, respectively. Retain 
only those arcs (u,w> for which (1) a + s + t < b , where s 

U U UV V u 

is the service time at customer u, and either (2) d + d < D, u,w 

U V 

e P , u * w, or (2’) d + d < D, u,w e P , u * w. Finally, 

u — n v - n 

let K be the fixed cost associated with each vehicle, incurred if 
this vehicle is utilized. 

Three types of variables are used in the mathematical 


formulation : binary flow variables , u , w N , v e V , u ^ w, 

time variables T , T 1 , T v , u e P, v e V, and load variables Y , 

u o 2n+i 

u € P. Let the binary decision variable equal to one if 

vehicle v travels from node u to node w, and equal to zero 
otherwise, v <s V, and u, w <= N. u * w. Let Tu be the time at 
which service at customer u begins, u € P. Let also be the 

time at which vehicle v leaves the depot and T^ n+1 be the time at 

which it returns to the depot, v € V. Next, let Y^ be the total 
load on the vehicle just after it leaves node u. u € P. It is 

assumed that the vehicles depart empty from the depot. i.e.. 


0 


constraints . 


These are decision variables with nonnegativity 


Finally, let f <T ) denote the penalty associated with 

beginning service at node u at time T^ , u € P, and g < zn+± o 

, , ~ r-rmtp of the indicated duration 

be the penalty associated with a rou 

. , . * i _ w *= u These penalties must be 

when traversed by vehicle v ^ 

. i. rnq fc r . We are now in a 

expressed in the same units as uw 
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position to present the pickup and delivery problem with time 
windows formally. The mathematical formulation is : 


Min £v<= vJ>=n £vew fu<Tu) + £vev g V (T^ n+i ~ T 0 > U> 






X v = 1 => Y - d = Y 

UW U V v 


> 

X 

= 1 *+ Y + 

d 

O v 

o 

V 

Y 

o 

= 0, 0 < Y 

< D 

x v 

UW 

b l nary , 



u e P , w <= P , 
v e V (1>> 

w e P + , v € V (16) 

U <= P + ( 17 ) 

u , w <= N , v <= V (18) 


We seek to minimize the sum of the total cost, the total 
penalty associated with servicing customers too early or too late, 
and the total penalty associated with routes exceeding a given 
duration. Constraints (2>-0> and (18) forma multicommodity 
min-cost flow problem. Constraints (6) ensure that the same 
vehicle, v, visits both u and n+u. Constraints (7) are precedence 
constraints which force node u to be visited befor node n+u. 
Next, constraints <8>-(10) describe the compatibility requirements 
between routes and schedules, while constraints (11)-(13) are the 
time window constraints. Finally, constraints (14)-(16) express 
the compatibility requirements between routes and vehicle loads, 
while constraints (17) are the capacity constraints. 


Note that the formulation includes a route duration 
restriction (b - a ). Note also that the formulation does not 

2rt+l O 

include subtour elimination. 

Let us now examin how models for different problem 
variants can be obtained from this formulation. The above 
forulat ion is the well known mu ltivehicle, man y to many 
dial-a-ride problem with time windows. If we let * V I = 1 and 
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e 1 imi nate constraint (6) we obtain the single vehicle variant. 

Redefine now N = P U CO.n+l), P = P + = p"= Cl,2,..,n}. 

To obtain the VRPTW we delete constraints (6>, (7) and (15). The 
constraints (8), as well as (14) ensure subtour elimination. If 

the VRPTW involves only pickup, define d > 0. For delivery 

V 

VRPTW, d < 0 , and in constraint (17), set Yo = D and eliminate 

V 

Y < D. The multiperiod VRP can be viewed as VRPTW where the 

U 

customers are replaced with request for service and for each 
request there is a time window. 

Furthermore, different schortest path problems with time 
window can be obtained from the above formulation. For example, 
constraints (1), (>)-<5), (8) — (15) , (18) and IVI = 1 constitue one 

shortest path problems with capacity, precedence and time window, 
can be obtained similarly. 
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PROGRAM FOR RAILWAY WAGON SCHEDULING 


PROGRAM heurilO (INPUT, OUTPUT); 

C Trying to determine costs. Used for solving problems 
with restricted capacity of arcs, so as to make mincost lnfil 
LABEL 25 , 50 , 75 , 100 , 200 , 250 , 500 , 550 , >75 , 580 , 550 , 575 , 580 , 600 ; 

CONST mmax=100; nmax=200; nodmax=75: arcmax=150; 

T_hormax = 120; C Max Time horizon (unit time) > 

TYPE 


C mainly used for the proce simplex ) 
matrix = ARRAY [l..mmax] OF REAL; 
ma t r i n t = ARRAY Cl..mmax] OF INTEGER; 
row__matrix = ARRAY Cl..nmax3 OF REAL; 
opt imi sat l ontype = ( maximisation, minimisation); 
bigmatrix = ARRAY Cl . .mmax , 1 . . nmax] OF REAL; 
inequalarray = ARRAY [l..mmax] OF CHAR; 

£ mainly used for defining node’s data } 
nodetype = 0..nodmax; 
nodeset = SET OF nodetype; 
nodrec = RECORD 

dema_supl : ARRAY [ 1 . . T_hormax] OF REAL; 
cum_dem_sup_g : REAL; 

requ_avai : ARRAY C 0 . . T_hormax3 OF REAL; 
cum_dem_sup_w : REAL; 

END; 

nodaray = ARRAY [l..nodmax] OF nodrec ; 

C mainly used for the defining arc's data } 
arc_rec = RECORD 

nodel , node2 : nodetype; 
d l st , cost_g , cost_w, 
capa, tmp_capa : REAL; 
de leted , enter : BOOLEAN; 

res_cap , 1 oad : ARRAY Cl - . T_hormax3 OF REAL; 
1 1 m s t , 1 1 m dt : REAL; 

END; 

arcaray = ARRAY Cl..arcmax] OF arc_rec; 



£ used for procedure f i nd_ass 1 gnment } 
dem_nod = RECORD 

nd t node type j 
dem r REAL; 
t lm ; 1 . . T__hormax ; 

END; 

ass_dem_nd = ARRAY Cl.. nodmax] OF dem_nod; 
supjiod a RECORD 

nd : nodetype; 
sup : REAL; 
t im ; 1 . . T__hormax ; 

END; 

ass_sup_nd = ARRAY [1 . .nodmax] OF sup^nod; 
sl_ass_typ = ARRAY Cl . . nodmax , 1 . . nodmax ] OF REAL; 
arcset = SET OF INTEGER; 
arcarr = ARRAY Cl.. arcmax ] OF INTEGER; 
del_type = ARRAY CO . . T_hormax] OF INTEGER; 
wtd_de lay_typ = ARRAY C 0 . . T__hormax ] OF REAL; 

VAR i mainly used for simplex procedure > 
mm, mml ,nnl ,ru , ii , j ) ,kk : INTEGER; 

Cm : row^ma trix; 

Am : b l gma trix; 

inequal itym ; i nequa 1 array ; 

bmi , bmo : matrix; 

ans : CHAR; 

p©riod_cost : real; 

inp.out ,out5 : TEXT; 

C u$ed in mincost network ) 

tt,T_, peri od_l an th , T_ hor : O..TJiormax; 
peri od_no s 1 . . 1 0 ; 

no^nod , no_arc, nod , nodi , nod2 : nodetype; 
arc : arcaray; 
node ; nodaray; 

degree; ARRAY Cl.. nodmax] OF INTEGER; 
ds : REAL; 

m i ncos t_f eas i b 1 e : BOOLEAN; 

C Note that since all the arcs undirected i.e. two way 
indegree = outdegree = degree /2 > 
cumsum : INTEGER; 

point, fmovptr , rmovptr : ARRAY Cl.. nodmax 3 OF INTEGER; 



f trace , rtrace : ARRAY Cl..arcmax3 OF INTEGER; 

C used for procedure f 1 nd_ass 1 gnment } 

dnn.snn, ino.jno, s_nod . cur_nod , d_nod : nodetype; 

dmn : ass_dem_nd ; spn : ass_sup nd; 

slak, assign : sl_ass_typ; 

dem_nod_set : nodeset; 

arc_flo_set : arcset; 

cn , ft , arc_s tor : INTEGER; 

assl , ass2 , ass , imn_arc_capa,prev_min_arc_cap, flo_sent, 
aval l_t im.mot i on_t ira, travl_t 1 m, s lack_t lm, de lay , rk REAL; 
s_tim,d_tim, time O..T_hormax; 

ddeg.sdeg. spt.dpt : ARRAY Cl..nodmax3 OF INTEGER; 
arcs_aray : arcarr ; 

C Empty wagon scheduling } 

no, i store , j store , kstore , ptr : INTEGER; 

excl_arc : arcset; 

maxload : REAL; passed, tried : BOOLEAN; 

C output information > 
delay_freq : del_type; 
wtd_de 1 ay : wtd_de 1 ay_typ ; 
tot_de lay , totai_cost : REAL; 
outl , out2 , out5 , out4 : TEXT; 

FUNCTION smaller <a,b ; REAL) : REAL; 

BEGIN 

IF a < b 

THEN smaller : = a 
ELSE smaller ; = b; 

END; 

PROCEDURE simplex 

t— 

<m,n: INTEGER; C: row_matrix; A:bigmatrix; 1 nequal l ty : l nequa larray ; 
VAR b_o: matrix; VAR dec i_vars : arcset ; VAR feasibility : BOOLEAN); 

LABEL 10,20* 

TYPE s e t_o f_var = SET OF 1 . . nmax ; 

VAR optimisation : opt 1 mi sat 1 on type ; 

set_of__s lack_var , set_of_art i f i c ial_var . set_of _a 1 trn_o P tm_ 

se t_of var ; 


b: matrix 


soln i 



degeneracy_f ound : BOOLEAN; 

PROCEDURE Resolve degeneracy ; 

c j 

VAR column, i2 : INTEGER; 

tempratio, rmnratio : REAL; 
degeneracy_resol ved : BOOLEAN; 

BEGIN C Resolve degeneracy > 

degene racy_re so 1 ved := FALSE; column := 0; 
WHILE ( NOT degeneracy_resol ved ) AND ( column <= m > DO 
BEGIN 

column : = column + 1; 
minratio := 10000; C initialise J 
FOR i2 := 1 TO stackptr DO 
BEGIN 

IF A_enter l ng[ l ndexroC l 2 3 3 > 0 

THEN tempratio := B_i nv[ l ndexroC i 2] , column ] / 

A_enteringC indexroC i233 

ELSE IF A_enteringC indexroC i233 = 0 
THEN tempratio := 10000 
ELSE tempratio := -10000; 

IF minratio > tempratio 
THEN BEGIN 

minratio := tempratio; 
indexrow : = indexro[i21; 

END; 

END; 

degeneracy resolved := TRUE; 

END C while } 

END; C Resolve degeneracy } 

BEGIN C Find min_ratio > 

degene racy_f ound = FALSE; 
min : = -1 ; 

FOR i 1 ; = 1 TO m DO 
IF rat i o[ i 1 3 >= 0 

THEN IF (min < 0) OR (min > ratio[il3> 

THEN BEGIN 

mm := ratioCil]; 
i ndexrow : = i 1 J 

degene racy__ found := FALSE; 



stackptr : = 1 . 



indexroCstackptr ] : = ,i ; 

END 

ELSE IF mm = ratioCil] 

THEN BEGIN 

degeneracy_found := TRUE; 
stackptr := stackptr + 1; 
indexroCstackptr] •= il- 

END; 

IF degeneracy_f ound = TRUE 
THEN Reso 1 ve_degeneracy ; 

END; C Find min ratio } 

PROCEDURE Int erchange_suf f 1 xes_of_var lab! es ( suff lxl , suf f 1 x2 : INTEGER); 
VAR temp ; INTEGER; 

BEGIN 

temp : = non_bas_var Csuf f 1 xl ] ; 

IF suffix2 <= m 
THEN BEGIN 

non_bas_var [ suf f i xl 3 : = bas l c_varC suf f l x2 ] ; 
bas l c_var[suf f l x2] := temp 

END 

ELSE BEGIN 

non_bas_varCsuf f i xl ] := non_bas_varCsuf f ix2] ; 

non_bas_var Csuf f i x2 ] : = temp; 

END; 

END; 

PROCEDURE Interchange (VAR x,y ; REAL); 

VAR temp ; REAL; 

BEGIN temp : = x; 
x : = y ; 
y ;= temp 

END; 

PROCEDURE I nterchange_co 1 umns_of_var tables 
(enter ing_bas i s , 1 eav l ng_bas l s : INTEGER); 

VAR il . INTEGER; 



VAR temp : REAL: 
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BEGIN 

FOR 1 1 := 1 TO m DO 
BEGIN 

temp := A[ i I , enter i ng_bas l s] ; 

ACil.enteri ng_bas l s ] : = A[ i 1 , 1 eav i ng_bas is]; 
AC i I , 1 eav i ng_bas l s ] := temp 

END 

END; 


BEGIN C Procedure simplex } 

REWRITE < out , ' outpro j e . dat ' ) ; 
pr i nt_no : = 1 ; 
ml : = m + 1 ; 

C To find out total no variables <nt) > 
nt := m+n; 

FOR l : = 1 TO m DO 
IF l nequal i ty C i ] = ’>' 

THEN nt := nt + 1; 
new_n : = nt ; 

optimisation : = minimisation; 

C Columns indexed m+1 to new_n are for non_bas_var 

and columns indexed 1 to m are basic_var } 

C Initialise non_bas__var and basic_var 1 

• FOR i ; B 1 TO nt-m DO non_bas_varCm+ 1 1 := l; 

FOR i s = 1 TO m DO basic_var[il ;= nt-m+i ; 

CASE optimisation OF 

maximisation i mi n_reducedcos t :=-l ; 

minimisation : max_reducedcost := 1; 

END; 

nl ; = m + n ; 

set_of_slac k_v a r : = CD; set_of_art i f i c l al_var ;= LI: 
S0<: _ o f_ a l trn_optm_soln s= C3; 

C Read fit 8e b } 

count := 0; 

FOR l ;=1 TO m DO 
BEGIN 

IF <nl+l <= nt) 
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THEN FOR j : = nl+1 TO nt DO 
AC 1 , ) 3 0; 

CASE mequal i tyC i ] OF 
f <’ i BEGIN 


+ Cbas i c_v arC i 3 3 


+ Cnon_bas_varCnl+count33 ; 


+ Cbas i c varC l 3 3 


C[ l 3 : = 0 ; 

s e t__o f _s 1 ack^var := set of slack var 


END; 

’>' : BEGIN 

count := count + 1 ; 

AC i , nl tcount 3 := -1; 

CCnl+count3 s - 0; 

set_of_s I ack__var : = set_of_s Iack_vai 

CASE optimisation OF 

maximisation : C[i3 := OM 3 -1000; 
minimisation : CC i 3 : = l M } 1000 
END; 

set of artificial var : = set of art 3 


END,; 

•= # : BEGIN 

CASE optimisation OF 

maximisation : CC l 3 s= C -M 1 -1 001 
minimisation : CC l 3 := C M > 1001 


END; 

set_of__art l f icial_var s = 

set_of_ar t 1 f i c i al^var *t Cbas i c_var C i J 3 

END; 

END; C case inequal i tyC i 3 > 

CBC i 3 := CC i 3 ; 

ENO; 


C Initialise coefficients of basic variables & B_mverse as unit matrix 

FOR l : = 1 TO m DO 
FOR j « a 1 TO m DO 
IF i - j 
THEN BEGIN 

AC i . j 3 s® 1 : 

B invC i . j 3 : = 1 
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END 

ELSE BEGIN 

At 1 , J 3 := 0; 

B_invC l , j 3 := 0 

ENO; 

C Initially find the SIMPLEX Multipliers } 

FOR j := 1 TO m DO 
BEGIN 

PIC jD := 0; 

FOR l :=1 TO m DO 

PICj] := PIC)] + CB[ i ] * B_inv[i,j] 

END: 

C Initially Find Reducedcosts for all variables 3 
FOR J : = 1 TO new_n DO 
BEGIN 

product := 0; 

FOR 1 : = 1 TO m DO 

product := product + PICi] * AC i , j 3 ; 
reducedcost C j ] ■ - Ct j ] - product 


ENO: 

C Initially take b_bar = b > 

FOR l := 1 TO m DO 
b_barCi] *= bCi]; 
i terat i on_no : - 0; deleted 0: 

CASE optimisation OF 

maximisation : F i nd_max_reducedcost < max_reducedcost 


pi vot__coi > : 


mi n imi sat i on 


Find min_reducedcost< mi n_reducedcost 


pivot_col ) ; 

END: 


5: REPEAT 

i terat ion_.no := i terat ion_no + 1 ; 

C Compute the updated coefficients of entering column ) 

FOR i s = 1 TO m DO 
BEGIN 

A_enteringC i ] = = 

FOR j := 1 TO m DO ■ 

A_enter i ngt i 3 == A_enter i ngC i 3 + B_mvCi.j3 

* AC j ,pivot_col 3 


END; 



o 0.001 


C Apply minimum ratio rule for deteriming leaving variable 
F i nd rat i os } 

TOR l s =1 TO m DO 
BEGIN 

IF abs < A_enter i ng[ i 3 - round (A_en ter mgC l 3 > > 


THEN A_enter i ng[ i ] := round(A_enter i ngC i 3 > ; 

IF A__enter i ng[ i 3 > 0 

THEN ratioCi] ; = b_bar[i3 / A_enter 1 ngC 1 3 
ELSE IF A_enter mgt 1 3 = 0 

THEN ratioCi] := 10000 
ELSE ratioCi] : = -10000; 

END; 

Find_min_ratio( min_ratio, pivot_row ); 

IF min_ratio < 0 
THEN BEGIN 

WRITELN ( out , 'This problem has 


unbounded solution ’>; 


WRITE ('Unbounded solution I!, 


Program aborted . - ' > ; 

GOTO 20 


C Update B 
FOR 


END; 

inverse matrix ! 

j . a l TO m DO CUpdate pivot row first! 

IF A_enter ingCpi vot_row] > 0 

THEN B mvCpi vot_row, ) 3 := B_t nvCp i vot_row . | 3 / 


A_enter i ngCp i vot__row] ; 

FOR i s = 1 TO m DO 


{ Update all other rows ! 


IF i <> pivot_row 
THEN FOR j:=l TO m DO 

B invC i , j 3 s = B_invC i , ) 3 


- 6 invCpi vot_row, ) 3 


* A_enter ingC i 3 ; 

{ Find b_bar ! 

FOR i t =1 TO m DO 
BEGIN 

b__bar C l 3 i = 0 • 

FOR J : = 1 TO m DO 
b bar[i3 : = b_bar[i3 


+ B l nvC l . J 3 * b[j3 


END; 

basic col i- pivot_row; 



£ Find S IMPLEX multipliers 'Pi's > 

CBCbas 1 c_co 1 3 : = C[p 1 vot_co 1 ] ; 

f OR j : a 1 ro m DO 

BEGIN 

PIC)] := 0; 

FOR i :=1 TO m DO 

PICj] := PIC,] + CBCi] * B_i nvC l , j ] 

END ; 

C Find Reducedcosts for all variables, 

note that reducedcost will automatically be = 0 

for bas l c_var > 

FOR , := 1 TO new_n DO 

BEGIN 

product := 0; 

FOR l : = 1 TO m DO 

product := product + PICiJ * ACi.j]; 
reducedcost C , ] := CC j ] - product 

END; 

Int erchange_suf f l xes_of_var l abl es (p i vot_coi , bas l c_col ) ; 
Interchange<CCpi vot_col ] , CCbas i c_col ] ) ; 

I nterchange_co 1 umns_of_var tables ( pi vot_col , bas i c_col > ; 
Interchange < reducedcost [pi vot_co 1 ] , reducedcost C bas i c_col 3 > 
pc := non_bas_var[pi vot_col 3 ; 

IF non_bas_var[p l vot_co 1 ] IN se t_o f_art if i c i a l_var 
THEN BEGIN 

£ Delete the artificial variable from subsequent considerations 
once it has left the basis } 

o 1 d_n : = nt - deleted; 

Int erchange_suf f i xes_o f_va nables(pi vo t__co 1 , o 1 d_n ) ; 
Interchange(C[pi vot_col 3 ,CCold_n] > ; 

Interchange_columns_of_var iab!es(pivot_col , old_n> ; 

Int erchange (reducedcostCpi vo t co 13# reducedcost Col d_n 3 ) ; 

deleted i- deleted + 1; 

END; 

feasibility := TRUE; 

FOR l : = 1 TO m DO IF b_barCi] < 0 
THEN BEGIN 

WRITELN (out, 'Now problem is primal infeasible, 

because one rhs constant b_bar is -ve’>; 

CASE optimisation OF 



/ 1 

maximisation : FOR j : = ml TO new_n DO 

IF not(j'IN set_of_s 1 ack_var ) 

AND not ( j IN a a t_o f _«i rtificial_var) and 

< reducedcost [ j ] > 0> 

THEN BEGIN 

WRITELN 

(out, 'It is also dual infeasible. Hence can ' ' t be solved further'); 

WRITE 

('Infeasible solution !!, program aborted 

feasibility := FALSE; 
GOTO 20; 

END; 

minimisation : FOR j : = ml TO new_n DO 

IF not(j IN set_of_slack_var) 

AND not ( ) IN se t_of_art 1 f 1 c i a l_var > and 

( reducedcos t [ j 3 < 0) 

THEN BEGIN 

WRITELN 

(out, 'It is also dual infeasible. Hence can ' ’ t be solved further’); 

WRITE 

('Infeasibible solution II, program aborted . . ’ ) ; 

feasibility := FALSE; 
GOTO 20 
END; 


END; 

WRITELN (out , 'However you may solve it further, 
by dual simlex method'); 

WRITE ('Primal infeasible !!, program aborted ..') 
GOTO 20; 


END; 

new_n ; = nt - deleted; 
CASE optimisation OF 

maximisation : BEGIN 


(max__reducedcost , pi vot_col ) ; 
(max__reducadcost> ) <= 0.001 
round(max reducedcost); 


F ind_max_reducedcost 
IF abs( max_reducedcost 
THEN max_reducedcost := 
IF max reducedcost = 0 


round 



THEN set_of_al trn optm_soln := 

set_of_altrn_optm_soln + Cpc] C Lnon_bas_v ar [ P i vot_col ] ] } ; 

END; 

minimisation : BEGIN 


(rm n_r educe dcost , p i vot_co 1 ) ; 
rouncKmi n__reducedcost ) > <= 0.001 
rouncKmin reducedcost ) ; 


E i nd_mi n_reducedcost 
IF abs ( mi n_reducedcost - 
THEN min reducedcost 


IF mi n_reducedcos t = 0 

THEN set_of_altrn_optm_soln : 

= set_of_aI trn_optm_soln + Cpc] C Cnon_bas_varCp i vot_co 1 ] ] ) ; 

END 

END; 

UNTIL < max__reducedcost <= 0 ) OR ( mi n_reducedcost >= 0 ) ; 

C Further minimisation of objective function is not possible 
i.e. when optimum solution is obtained }; 

IF set_of_al trn_optm_soln <> [] 

THEN BEGIN 

WRITE <out,'Set of alternate optimum soln contains : '); 

FOR i : = 1 TO nt DO 

IF l IN set_of_al trn_optm_soln 

THEN WRITE ( ou t . ' X ' , I : 2 ) ; WRITELN (out); 

END; 

10: C Compute the value of minimum value of objective function} 

Z := 0; 

, FOR i : = 1 TO m 00 

Z : = Z + CBCil * b_bar [ l ] ; 

C Print output > 

CASE optimisation OF 

maximisation : WRITELN (out , 'Maximum objective function value 
is Z = * , Z : 6 : 2 ) ; 

minimisation : WRITELNfout . ‘Minimum objective function value 
is Z = ’ ,Z;6:2> 

END; 

FOR i := 1 TO m DO 

b_o C i 3 := 0; C initialise before returning output > 

WRITELN(out , ’Opt imum solution is ’>: 

FOR i := 1 TO m 00 



BEGIN 


IF NOT ( <bas i c_var[ 1 ] IN set_of_s lack_var ) OR 
(bas 1 c_varC i 3 IN set_of_art i f i c ial_var ) ) AND <b_bar[i] > 0.001) 

THEN BEGIN 

WRITE (out.* X ’ , bas l c_var [ l 3 : 2 , ' = * . b_barC l 3 : 4 : 2 > ; 
deci_vars := deci_vars + Cbas i c_var[ l ] 3 ; 
b_ot bas i c_var[ i ] ] := b_bar[i3; 

END; 

END; WRI TELN< out ) ; 

FOR i := 1 TO m 00 

IF ( bas i c_var C l 3 IN se t_of_art l f i c l al_var ) AND (b_barCi3 <> 0 ) 
THEN BEGIN 

WRITELN < out , 'Since artificial variable is present 

in the optimum solution '); 

WRITELN (out . 'PROBLEM IS INFISIBLE'); 

WRIT ELN(out , ' X* ,basic_var[ l 3 :2 , ' = ’ . b_bar[ i 3 : 4 : 2 > ; 

WRITELN ( ' Inf i s lbi 1 l ty !!. program aborted ..'); 
feasibility :== FALSE; 

GOTO 20; 

END 

ELSE feasibility := TRUE; 

CASE optimisation OF 

maximisation : IF (max_reducedcost = 0) 

AND not(non_bas_varCpivot_col 3 IN set_of_altrn_optm_soln> 

THEN BEGIN WRITELN (out. 

•Finding out another alternate optimum soln ’>;GOTO 5 END; 

minimisation : IF (mi n_reducedcost = 0) 

AND not (non_bas_varCpi vot_col 3 IN set_of_al t rn_optm_so 1 n ) 

THEN BEGIN WRITELN ( out . ’ F l nd i ng out 

another alternate optimum soln ');G0T0 5 END, 

END; 

WRITE ('Complete solution has been written in the 

file outproje.dat'); 

20: READLN ; (CLOSE ( out >; > 

END; C Procedure Simplex > 

FUNCTION shortest_dist ( s. t : nodetype): REAL; ^ 



VAR maxdist, temp_dist ; REAL; 

Neet, temp, perm : nodeset; 
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shordist : ARRAY Cl..nodmax] OF REAL; 
i.j, arc_no : INTEGER; 
store, storl : nodetype; 

FUNCTION min_temp_dist_nod : nodetype ; 

VAR minimum : REAL; k : nodetype; 

BEGIN 

minimum : = maxdist; 

FOR k := 1 TO no_nod DO 
IF NOT <k IN perm) 

THEN IF shordistCk] < minimum 
THEN BEGIN 

minimum := shordistCk]; 
min_temp_di st_nod := k; 

END; 

END; C mi n_temp_di st_nod > 

BEGIN C function shortest dist } 

Nset := CI..no_nod]; 

perm := Cs3; temp := Nset - Cs ] ; 

IF NOT (<s IN Nset) AND <t IN Nset)) 

THEN BEGIN 

WRITELN < ’Node out of network data range !!’>; 

GOTO 500; 

END; 

£ Set shortest dist for other nodes = sum of all dist.s in network > 
maxdist := 0; 

FOR i : = 1 TO 2*no_arc DO 

maxdist := maxdist + arcCil.dist; 

FOR i := 1 TO no_nod DO 

shordistCi] := maxdist; 
shordistCs] s= 0; 
store := s ; 

WHILE NOT <t IN perm) 00 
BEGIN 

FOR j : = pointCstore] TO po i nt [store+1 3 -1 DO 

BEGIN C Get the node in the adjuscency list of node -store- 
called * storl * , update its shordist if required > 
arc_no := ftrace[j3; 

t emp_d l s.t ;= shordistCstore3 + arcCarc_no3 . di st ; 



storl := arc[arc_no] .node2 ; 

IF storl IN temp 

1MFN shordistCstorl ] t « smaller < shord i ® t C ® tor 1 ] * 

temped l st ) ; 

END; C for adjuscent nodes update shordist ) 
store := mi n_temp_di st_nod ; 
temp := temp - [store]; 
perm := perm + [store]; 

END; C while not (t in perm) DO } 
shortest_di st := shordist[t]; 

END; C function Shortest dist ) 

PROCEDURE Find_Assignment < nodn.nosn : nodetype; 

^ dn : ass_dem_nd; sn : ass_sup_nd; 

slack : sl_ass_typ; 

VAR assignment : sl_ass_typ); 

VAR ip. jp. ro.co, count ; INTEGER; 

assigned : ARRAY [1 . . nodmax , 1 . . no'dmax] OF BOOLEAN; 
totslack, small : REAL; 


PROCEDURE Max_slack <VAR row, col: INTEGER); 

VAR i, j . INTEGER; max : REAL; 

PROCEDURE braketie; 

VAR penaltyl, pona 1 ty E : REAL; 

FUNCTION next_max_in_row (iro : INTEGER; : REAL. 

VAR k ; INTEGER; nextmax : REAL; 

BEGIN 

nextmax : = 0; 

FOR k := 1 TO nodn DO 

IF (nextmax < slackC iro.kj) AND CUcKC . ro . K] < -»»> 

AND NOT ass i gnedC i ro , k ] 

THEN nextmax : = s lackC i ro , k] ; 
next max_in_row : = nextmax; 


END; 

FUNCTION next_max_in_col <ico 
VAR k : INTEGER; nextmax 

BEGIN 


INTEGER) : REAL; 
REAL; 


nextmax : = 0; 


FOR k x= 1 TO nosn DO 
IF (nextmax < s lackCk . i co] > AND 


(sl ack[k 


co ] < max) 


/ 
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AND NOT assignedtk, ICO] 

THEN nextmax := slacktk , i co] ; 
next_max_tn_col : = nextmax; 

END; 

BEGIN C braketie } 

IF j = col 
THEN BEGIN 

penal tyl : = max - next_max_i n_row( 1 ) ; 
penal ty2 := max - next_max_i n_row( row) ; 

IF penal tyl > penalty2 
THEN row : = 1 ; 

END; 

IF i = row 
THEN BEGIN 

penal tyl := max - next_max_i n_co 1 ( j ) ; 
penalty2 : = max - next_max_i n_co 1 (col > ; 

IF penaltyl > penalty2 
THEN col := j ; 

END; 

END; C braketie } 

BEGIN C procedure max_slack } 

, max := 0; row := 1; col := 1; 

FOR i := 1 TO nosn DO 
FOR j := 1 TO nodn DO 

IF Cmax < slackCi.jl) AND (dn[j].dem> 0.001) AND (snCil.sup > 
THEN BEGIN 

max : = slackC l . j ] ; 
row : = i ; 
col : = j ; 

END 

ELSE IF ( round(max) = round! s lackC l . j 3 > > AND 

( dnC j ] . dem > 0.001) AND (sn[i].sup > 0.001) 

THEN IF NOT ass i gned[ l . j ] 

THEN IF < l =row) OR Cj=col) 

THEN braketie 
ELSE BEGIN 

row : = i ; 
col : = ) ; 

END; 

END; C procedure max_slack > 

/ 


0 . 001 ) 



FUNCTION demand_satisf led : BOOLEAN; 

VAR 1 i INTEGER; 

BEGIN 

domand_aat 1 sf 1 ed := TRUE; 

FOR l : = 1 TO nodn DO 
IF dn[ i ] . dem > 0 . 001 
THEN demand_sat l sf i ed := FALSE; 

END; 

PROCEDURE pr i nt_ass l gnment ; 

VAR l.J : INTEGER; 

BEGIN 

FOR l : = 1 TO nosn 00 
BEGIN 

FOR j : = 1 TO nodn DO 
WRITE (ass l gnment [ l ,j]: 4 : 0 ) ; 

WRI7ELN; 

END; 

READLN ; 

END; 

BEGIN C Procedure F l nd_ass l gnment } 

FOR ip : = 1 TO nosn DO 
FOR jp ;= 1 TO nodn DO 
BEGIN 

assignmentCip.jp] := 0; 
assignedCip.jp] := FALSE; 

END; 

count := 0; 

REPEAT 

count ;= count + 1; 

Max_s lack (ro.co); 

small := smai ler(snCro] .sup.dnCco] .dem > 

snCro] . sup : = snCro].sup - small; 
dnCco] .dem : = dnCco].dem - small; 

Ass l gnment C ro , co] := small; 
ass i gnedC ro , co] := TRUE; 

£ pr i nt_ass l gnment ; 1 

UNTIL <demand_sat i sf l ed) OR (count >= 200); 



can * # t be satisfied 


IF count = 200 THEN WRITELN (’Demand 
tots lack := 0 ; 

FOR ip : - ] 7 0 nosn DO 
FOR jp := I TO nodn DO 

tots lack := tots lack + slackC ip, j p]*ass i gnment C ip, jp] 
WRITELN ('Total slackness : tots lack :$■ 1 ) • 

END; C Procedure F l nd_ass l gnment > 


PROCEDURE Find_short_assign < nodn. nosn : nodetype; 

^ ^ * ass_dem_nd ; sn : ass sup nd; 

short : sl_ass_typ; 

VAR assignment : sl_ass_typ>; 

VAR ip. JP. ro.co. count : INTEGER; 

assigned : ARRAY Cl . . nodmax , 1 . . nodmax ] OF BOOLEAN; 
totshort, small : REAL; 

PROCEDURE Min_short (VAR row, col: INTEGER); 

VAR i.j : INTEGER; min : REAL; 

PROCEDURE braketie; 

VAR penaltyl, penalty2 : REAL; 

FUNCTION next_min_in_row ( i ro : INTEGER) : REAL; 

VAR k : INTEGER; nextmin : REAL; 

BEGIN 

nextmin := 1000; 

FOR k := 1 TO nodn 00 

IF (nextmin > shortC l ro , k] ) AND ( short C i ro , k ] > min) AND 
NOT ass i gnedC l ro , k] 

THEN nextmin := short C i ro . k 3 ; 
next_min_in_row := nextmin; 

END; 

FUNCTION next_mi n_i n_col < i co : INTEGER) : REAL; 

VAR k ; INTEGER; nextmin : REAL; 

BEGIN 

nextmin := 1000; 

FOR k := 1 TO nosn DO 

IF (nextmin > shor t Ck . l co] ) AND ( short Ck . l co] > mm) 

ANO NOT ass i gnedCk . l co] 


THEN nextmin := short C k . i co 3 : 
next min in col := nextmin; 


END; 

BEGIN C braketie > 


/9 


IF j = col 
THEN BEGIN 

penalty! := next_min_in_row< i ) - mm- 
penal ty2 ;= next_rm n_i n_row< row) - min; 

IF penal tyl > P enalty2 
THEN row : = i ; 

END; 

I F i = row 
THEN BEGIN 

penal tyl : = next_mi n_i n_col < , ) - min; 
penal ty2 := next_mi n_i n_col ( col ) - min; 

If- penal tyl > penalty2 
T HEN col : = j ; 

END; 

END; C braketie } 

BEGIN C procedure Min_short } 

min := 1000; row := 1 ; col := 1; 

FOR 1 : = 1 TO nosn DO 
FOR J ; * 1 TO nodn DO 

IF <min > shortCi.j]) AND < dn[ j ] . dem> 0.001) 

AND (sn[ i ] .sup > 0.001) C preferably check time feasibility 1 
THEN BEGIN 

min : = shortC l , j ] ; 
row .- = i ; 
col : = j ; 

ENO 

ELSE IF (round(min) = round! short [ i , j 3 ) ) AND 

(dn[j].dem > 0.001) AND (sn[i].sup > 0.001) 
THEN IF NOT ass i gnedC i , j] 

THEN IF <i=row) OR <j=col) 

THEN braketie 
ELSE BEGIN 

row : = i ; 
col : = J ; 

END; 

END; C procedure min_short 1 


FUNCTION demand satisfied : BOOLEAN; 



ou 

VAR l : INTEGER; 

BEGIN 

dcmand_satisf ied := TRUE; 

FOR 1 := I TO nodn DO 

IF dnC i ] . dem > 0.001 
THEN demand_sat isf ied := FALSE; 

END: 

PROCEDURE pr l nt_ass i gnment ; 

VAR l , j ; INTEGER; 

BEGIN 

FOR l := 1 TO nosn DO 
BEGIN 

FOR j ; = 1 TO nodn DO 
WRITE <ass l gnment [ i , j ] : 4 : 0 ) ; 

WRITELN; 

END; 

READLN ; 

END; 

BEGIN C Procedure F i nd_short_ass i gn 1 
FOR ip : = 1 TO nosn DO 
FOR jp := 1 TO nodn DO 
BEGIN 

assignmenttip.jp] ; = 0; 
assignedCip.jp] : = FALSE; 

END; 

count ;= 0; 

REPEAT 

count : = count + 1 ; 

Mi n_short (ro.co); 

small := smal ler(snCro] .sup.dnCco] .dem >: 
snC ro] . sup : = snCroD.sup - small; 
dnCco].dem := dnCcol.dem - small; 

Ass i gnment C ro.co] := small; 
ass i gnedC ro , co ] := TRUE; 

C pr l nt_ass i gnment ; } 

UNTIL <demand_satisf ied) OR (count >= 200); 

IF count = 200 THEN WRITELN < ’Demand can * ’ t be satisfied **’ > 


totshort := 0; 
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FOR ip : = 1 TO nosn DO 
FOR )p := 1 TO nodn DO 

totshort ,= totshort t shor t [ , p . )P 3*a SS , gnment C , p , , p] , 
WRITELN ('Total shortness : *, totshort : 5 : 1 ) : 

END; C F i nd_short_ass i gn > 

PROCEDURE f ind_short_ex_aray ( s. t : nodetype; ex_arc : arcset; 

VAR short_aray ; arcarr; VAR stackptr : INTEGER ); 

£ } 

LABEL 20; 

VAR maxdist, temp_dist : REAL; 

Nset . temp, perm : nodeset; 
shordist : ARRAY [ 1 . . nodmax] OF REAL; 
i,j, arc_no : INTEGER; 
store, storl : nodetype; 

FUNCTION mi n_temp_d i st_nod : node type ; 

VAR minimum : REAL; k : nodetype; 

BEGIN 

minimum := maxdist; 

FOR k i= 1 TO no_nod DO 
IF NOT <k IN perm) 

THEN IF shordistCk] < minimum 
THEN BEGIN 

minimum ;= shordistCk]; 
mi n_temp_d i st_nod := k; 

END; 

END; C mi n_t emp_d i s t_nod } 

BEGIN C Procedure find short aray } 

Nset := Cl.. no_nod] ; 

perm : = C a ] ; temp :* Nset - Cs]; 

IF NOT <<s IN Nset) AND (t IN Nset)) 

THEN BEGIN 

WRITELN ('Node out of network data range !!'); 

GOTO 500 ; 

END; 

C Set shortest d.st for other nodes = sum of all dist.s in network 
maxdist : = 0: 


FOR i:- 1 TO 2 *no_arc DO 



maxdist : = maxdist + arc[i].dist; 

FOR 1 : = 1 TO no_nod 00 

shordistti] ■= maxdist; 
shordistCs] := o ; 
store : = s ; 

WHILE NOT (t IN perm) DO 
BEGIN 

FOR J: = po i nt [ store ] TO po i nt [ store+1 3 -1 DO 
BEGIN C Get the node in the adjuscency list of node ’store' 
called storl ' , update its shordist if required } 
arc_no := f traced j] : 

IF not (arc_no IN ex_arc) 

THEN temp_dist : = shordistCstore] + arc[arc_no! . di s 
ELSE temp_dist := 10000; 
storl := arc[arc_no] . node2 ; 

If storl IN temp 

THEN shord l st[ storl ] := smaller (shordisttstorl 3 , 

t emp_d l s t > ; 

END; C For adjuscent nodes update shordist 3 
store := mi n_tamp_di s t_nod ; 
temp ;= temp - [store!: 
perm := perm + [store]; 

END; C while not (t in perm ) DO } 

C find the arcs in the shortest path } 

perm := C 1 3 ; temp := Nset - [t]; 

* 

store := t; stackptr : = 0; 

I WHILE NOT (s IN perm) DO 
BEGIN 

FOR )!* po i nt [ store 3 TO po i nt [ store + 1 3 -1 DO 
BEGIN 

arc_no : = rtraced j 3 ; 

f IF < shord i st [ arc [arc_no3 . nodel 3 = shord l st[ store! 

- arc[arc_no 3 . d i st ) AND NOT (arc_no IN ex_arc> 

THEN BEGIN 

store := arc[arc_no! .nodel ; 

GOTO 20; 

END; 

END; 


20 i 


stackptr i * stackptr + 1.* 
short arayCstackptr] := arc_no; 



temp := temp - [store]; 
perm := perm + [store]; 

END; C while not (s in perm ) do } 

END; £ Procedure find short ex aray } 

PROCEDURE Update_de lay__i nf o(de 1 e : REAL; flo : REAL; 

£ 

VAR de l_aray : del_type; VAR wt_dele : wtd delay typ; 

VAR stackptr : INTEGER); 

VAR 1 nt_de 1 : INTEGER; 

BEGIN 

int_del := round(dele); 

IF int_del > stackptr 
THEN stackptr := int_del; 

de l_aray C i nt_de 1 ] := de l_aray C i nt_de 1 ] + round(flo); 
wt_de 1 e [ l nt_de 1 ] : = wt_de 1 e [ i nt_de 1 ] + dele*flo; 

END; 


BEGIN C main i 

WRITE (’Do you want to begin with infinite capacity of arcs? ') 
READLN Cans); 

IF (ans = ' y ' ) OR (ans = 'Y' ) 

THEN RESET £ i np , ' i npsup . dat ’ ) 

ELSE RESET < inp, * mpsupcap.d' ) ; 

READLN ( i np , no_nod , no_arc); 

FOR i i : = 1 TO no_arc DO 

WITH arcC i i ] DO 

IF (ans = ' y ' > OR (ans = ’Y'> 

THEN READLN (inp, nodel , node2.dist) 

ELSE READLN (inp. nodel, node2.dist.capa); 

FOR n : * 1 TO no_arc DO 

WITH arcC i i 3 DO 

BEGIN 

cos t_g := dist / 10; 

cost w ; = cost_g /l. 5; 

dist := dist / 00*4); 

IF (ans = ’ y ' > OR (ans = ’Y’) 


/ 


THEN capa •• = 1000; 
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ENO; 

FOR 11 := no_arc+l TO 2*no_arc DO 

BEGIN 

arcC 1 l 3 . nodel : = arc[ i i -no_arc] . node2 ; 
arcC i i 3 . node2 := arc[ l i - no_arc] . nodel ; 
arcC i i 3 . cost_g := arcC l l -no_arc 3 . cost_g ; 
arcC l i 3 . cost_w := arcC l l -no_arc3 . cost_w; 
arcCii3.dist := arcC l i - no_arc3 . dist ; 
a rc C i i 3 . capa : = arcC i i - no_arc3 . capa; 

ENO; 

WRITE (’Give period length ; ’ ) ; READLN (per i od_l enth) ; 
WRITE (’Give time horrizon : ’>; READLN (T_hor>; 

IF T__hor < period_lenth 
THEN BEGIN 

WRITELN (’Time horrizon must be > period lenth ’> 
GOTO 25; 

END; 

FOR ii := 1 TO 2*no_arc DO 
WITH arcC i i 3 DO C initialise 1 
FOR tt := 1 TO T_hormax DO 
BEGIN 

r es_capC 1 1 3 : = capa; 
loadCtt3 ;= 0; 

ENO; 

FOR ii : = 1 TO no_nod DO 
WITH nodeCii] DO £ Initialise 1 
FOR tt ;= 1 TO T_hor DO 
BEGIN 

dema_supl C 1 1 3 ; = 0; 

requ_avai [tt] ;= 0: 

ENO; 

WHILE NOT EOF (inp> DO 
BEGIN 

READ (inp. nod ,ds , time) ; 

WITH nodeCnod] DO 


dema_supl C t ime*6 3 == ds 

IF EOLN ( mp) 

THEN READLN(mp): 


ENO; 
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RESET ( i np , ' inpdem. dat ' > ; 

WHILE NOT EOF (mp) DO 
BEGIN 

READ (mp, nod, ds, time) j 
WITH nodeCnod] DO 
dema_supl C t ime*6 3 : = -ds ; 

IF EOLN (inp) 

THEN READLN(inp); 

END; 

RESET ( mp. ’ mprowag.dat’ ) : t Check and delet to read i 

WHILE NOT EOF (inp) DO 

BEGIN 

READLN < i np, nod, ds , t ime) ; 

WITH nodeCnod] DO 

requ_avai Ctime*6] := requ_avai [t ime*63 + ds ; 

ENO; 

CLOSE (mp); 

FOR i l : = 1 TO no_nod DO 
WITH nodeC i i 3 DO 

cum_dem_sup_g i= 0; 

FOR li : = 1 TO no_nod DO 
degree C i l 3 : = 0 ; 

FOR ii : = I TO 2*no_arc DO 
WITH arcC i l 3 DO 

degreeCnode23 := degree[node2] + 1; 

cumsum : = 1 ; 

FOR n != 1 TO no__nod+l 00 
BEGIN 

pomtC l i 3 cumsum; 

cumsum := cumsum + degreeCii]; 

ENO; 

FOR 11 ; » I TO no_nod+l DO (Initialisation > 

BEGIN 

fmovptrCii] := pointCii]; 
rmovptrCii] ;= pointCiiJ; 


END; 

FOR ii ;= 1 TO 2*no_arc DO 

WITH arcC 1 1 3 DO 

BEGIN 


t once > 


/ 



ftraceC fmovptrlnodel ] ] .= (I . 
fmovptr [nodel 3 : = fmovptr [nodel ] + j. 
rtracaC rmovptr [node 2 ] ] ,= , j . 
rmovptr[node23 : = rmovptr Cnode2 ] + i, 

END; 

FOR 11 : = 0 TO T_hor DO 

BEGIN 

del ay__f req[ 1 i ] : = 0; 
wt d__de 1 ay [ 1 i ] : = Q; 

END; 

REWRI TE <out2, ' out route . dat * ) ; 

REWRITE < out 5 , 1 outcost . dat * ) ; 

WRITELN (out?, 'Period Goods Cost Empty Wagon Cost’); 
par l od_no i= 1 ; 

WRITELN; WRITELN ('Goods wagon scheduling'); 

WRITELN ('Period no : ' , per i od_no ; 2) ; 

T_ := period_no* per l od_l enth ; 

FOR i i : = 1 TO no_nod DO 
WITH nodeC l l 3 DO 
BEGIN 

FOR tt := (per iod_no-l )*per iod_lenth + 1 

to period_no*per iod_lenth DO 
cum_dem_sup_g : = cum_dem_sup_g + dema_supl [ tt 3 ; 
IF cum_dem_sup_g < -0.001 
THEN dem_nod_set : = dem_nod_set + Cii3 

END; 

mm no__nod + no_arc *2; C No. of constraints } 
ni s* 2*no_arc; C No of decision variables > 

FOR ii :* 1 TO 2*no_arc DO 

cmCmm+ii3 :* arc! i i 3 . cost_g; 

FOR n « * 1 TO mm DO 
FOR jj ;= mm+1 TO mm+m DO 
AmC i l , j j 3 ; = 0 ; 

FOR i i ; * 1 TO no_nod DO 
WITH node! l i 3 DO 
BEGIN 

C For all the arcs incoming at node ii > 

FOR jj ;= pointCii3 TO pointCn+13 -1 DO 


BEGIN 
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kk : = mm + rtraceCj]3; 

IF cum_dem__sup g <= 0 
(HEN AmC i i , kk ] : = 1 
ELSE AmCii.kk] •.= -l; 

END j 

C For ail the outgoing arcs from node n } 
FOR jj := pointCii] TO pointCii+1] -1 DO 
BEGIN 

kk := mm + f traceC j j 3 ; 

IF cum_dem_sup_g < = 0 
THEN AmC i i .kk] -1 
ELSE AmC l l , kk] : = 1 ; 


END,- 

IF cum_dem_sup_g < -0.001 
THEN i nequal l tymC i i 3 := ’>' 

ELSE IF abs(cum_dem_sup_g) <= 0.001 
THEN l nequal i tymC i i ] ;= ’ = ' 
ELSE inequal i tymC i i 3 := 
bmiCii3 := abs(nodeC l l 3 .cum_dem_sup_g) ; 
END; C for li ; = 1 TO no_nod do 1 
WRITELN (‘Node cum_dem_sup_g ' > ; 

FOR li ; = 1 TO no_nod DO 

WRITELN ( i l : 4, nodeC i i 3 . cum_dem_sup_g : 8 : 1 > ; 

Jj := 0; 

FOR n ;= 1 TO no_arc*2 DO 
BEGIN 


J j != 1 j + 1 J 
AmC i i +no_nod ,mm+ j j ] ••= 1; 

i nequal i tymC i i +no_nod] ;= 

bmi Cno_nod+ 1 l ] : 3 arcC i i 3 • capa ; 


END; 

simplex <mm,ni , Cm, Am. i nequal ltym.bmi 
period_cost s= 0; 

IF mi ncost_f eas lbl e THEN 

FOR ii := 1 TO 2*no_arc DO 
per l od_cost t= period_cost + 
snn ;= 0; C No. of sup nodes ) 
dnn ;= 0; C No. of dem nodes J 


bmo , arc_f lo_set ,mi ncost_f eas i bl e ) 


bmoC i i 3*arcC i i 3 - cost_g; 


FOR ii ; = 1 TO no nod+1 DO 


RPP IM 



ddegl i 1 j : = 0 ; 
sdegC 1 1 3 : = 0 ; 


END; 

FOR i i : = 1 TO no_nod DO 
WITH nodeC i i 3 00 
BEGIN 

FOR tt := I TO T_ DO 
IF dema_supl C 1 1 3 < -0.001 
THEN BEGIN 

dnn : = dnn + 1 ; 
ddegCii3 := ddeg[ii3 + 1; 

WITH dmn [dnn 3 DO 
BEGIN nd := n ; 

dam : = abs (dema_supl C tt 3 > : 
tim := tt; 

END; 

END 

ELSE IF dema_supl C tt 3 > 0.001 
THEN BEGIN 

snn : = snn + 1 ; 
sdeg[ii3 : = sdeg[i;3 + 1; 

WITH spn[snn3 DO 
BEGIN nd := l i ; 

sup ;= dema_supi Ctt3 ; 
t im : = tt ; 

ENO; 

END; 

ENO; 

IF mi ncost_f eas i bl e 
THEN BEGIN 

FOR i i s= 1 TO snn DO 
FOR J J : = 1 TO dnn DO 
BEGIN 

i no : = spnC i l 3 . nd ; 
j no : = dmnC j j 3 . nd ; 

s lakC i i . ) J 3 := dmn [ jj 3 . 1 1 m - spntii3.tim 

- short est_d l s t < l no , j no ) 


END; 

f i nd_ass i gnment ( dnn . snn . dmn . spn . s 1 ak , ass i gn ) 


ENO 
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ELSE BEGIN 

FOR it : = 1 TO snn DO 
FOR jj ;= 1 TO dnn DO 
BEGIN 

i no s = spnC i t ] . nd ; 
j no : = dmnC ) ) ] . nd : 

slak[it,jj] •= shor tea t_d t 3 t ( t no , j no) t 

END; 

F « nd_shor t_as s tgn<dnn, snn, dmn , spn , s lak , ass i gn> ; 
FOR li := I TO 2*no_arc DO 
WITH arc [ t t ] DO 

FOR tt := (pert od_no- 1 ) #pe r i od_l enth + 1 TO 
per t od_no#per t od_l enth DO 
res_cap[tt] := capa/6; 

ENO; 

FOR tt := 1 TO 2*no_arc DO 

WITH arcC i i ] DO 

BEGIN 

tmp_capa : = bmoC it]; 

IF tmp_capa < 0.001 
THEN deleted := TRUE 
El.SE deleted := FALSE; 

ENO: 


cumsum t~ 1 ; 

FOR ii := 1 TO no_nod+l DO 
BEGIN 

dpt C i i ] : = cumsum; 

cumsum := cumsum + ddegCii]; 

END: 

cumsum : = 1 : 

FOR 1 i := 1 TO no_nod +1 DO 
BEGIN 

sptt i t ] := cumsum: 

cumsum := cumsum + sdegtit]: 

ENO: 

FOR i i : = 1 TO dnn DO 

WRITE (dmnC i i3 . nd: J > t 



WRITELN; 


WRITE (’ '); 

FOR 11 : - 1 10 dnn DO 

WRITE <dmnC i 1 3 . dem: 5 : 0 ) j 
WRITELN; 

WRITE (’ •); 

FOR ii : = I TO dnn DO 

WRITE (dmnC l l 3 . t im: > ) ; 

WRITELN; 

FOR s_nod : = 1 TO no_nod DO 

FOR j j i = sptCs_nod3 TO sptCs_nod+l 3-1 DO 

BEGIN 

WRITE ( spnC ) ) 3 . nd : J , spn[ j j 3 . sup : J ,• 0 , spnC j j 3 
FOR cur__nod := 1 TO no_nod DO 
FOR kk : = dptCcur_nod3 TO dpt[cur_nod+l 3-1 DO 
WRITE (ass 1 gnC j j , kk3 : 5 : 0 ) ; 

WRITELN; 

END; 

WRITELN (out2, ’Goods Period no : ' , per 1 od_no : 2 > ; 

WRITELN (out 2 .’ Flow Node Time 

cri : 3 1 ; s_nod : = 1 ; 

IF mi ncost_f eas ibl e THEN 

WHILE (dem_nod_set <> C3> AND <cn<400> DO 
BEGIN 

WHILE <nodeCs_nod3 . cum_dem_sup_g < 0.001) 

AND (s_nod < 3 no_nod) DO 

s_nod := s_nod + 1 ; 
cur_nod ; = s_nod; 

WRITELN Ocn ; ',cn:2,' s_nod'. s_nod:4,* 

' , node C s_nod 3 . cum_dem_sup_g 1 4 : 0) ; 

min_arc_capa := 1000 ; 
no : = 0 ; 

t ra v 1 __t 1 m s = 0 ; flo_sent : = 0 ; 

FOR ii := 1 TO 2 *no_arc DO 

arcC i i 3 . enter : - TRUE ; 

7^.- REPEAT 

WRITELN <’cur_nod : ' , cur_nod : 4> ; 


t imi 5 


> ; 


sup 


WRITE (’ft 

FOR ii P ointCcur_nod3 TO po.ntCcur_nod+l 3-1 DO 



WITH arcC f traceC i i ] 3 DO 
BEGIN 




nun_arc_capa> ; 


200 : 

t im-spnC j store ] 

inf eas ibi 1 i ty 
- aval l_t i m; 
delay : 5 : Z ) ; 


ft := f t race [ i i ] ; WRITE (ft:4); 

IF (ft IN arc_f 1 o_se t ) AND NOT deleted AND enter 
THEN BEGIN 

WRITELN ( ' enters f ) ; 

arc_stor : = ft; 

no : “ no f 1 ; 

arcs_aray[no] : = ft; 

travl_tim : = travl__tim + dist; 

pre v_mi n_arc:_ ca P : - rm n_arc_capa ; 

mi n_arc__capa : = smal 1 er (arcC f 1 3 . tmp_capa , 

IF node2 IN dem_nod_set 
THEN BEGIN 

ass := 0; 

FOR j j : = spt [ s_nod] TO spt[s_nod+l 3~1 00 
FOR kk : - dpt[node23 TO dpt[node2+l 3-1 DO 
IF ass ign[ ) ) ,kk3 > 0.001 
THEN BEGIN 

ass := assignC ) j .kk3 ; 
j store : = j J i 
kstore : = kk ; 

GOTO 200 ; C ex 1 1 for loop> 

END; 

WRITELN ('Demand nod ' , node2 : 4) ; 

IF ass > 0.001 
THEN BEGIN 

ava i l_t 1 m : = dmnCkstore3. 

IF avail_tim < travl_tim 

THEN BEGIN 

WRITE ('time 

delay := travl_tim 

WRITELN ( 'delay ' , 

s tim := spnCjstore 


tim. 
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round< travi_t im> ; 


d tim : = s t 1 m + 


> 0.001 


END 

ELSE BEGIN 
cL.tim : = dmnCkstore] . t im; 

delay : = 0 ; 

IF 7 d_tim - round(travl tim) 


-round < t rav 1 t lm) 


THEN s tim := d tim 


de lay_f req , wtd_de lay ,ptr>; 
- f 1 o sent ; 


ELSE sjim := 0; 

END; 

flo_sent := sma 1 1 er (mi n_arc_capa , ass); 
Updat e_de lay_i nf o ( delay, f lo_sent , 

ass i gn[ j store , ks tore 3 := assignC}store,kstore] 

WITH dmnCkstore] DO 
BEGIN 

dem : = dem - fio_sent; 
FOR jj := 1 TO snn DO 
assign C J ) .kstore] : ~ 


ass ignC j J . kstore] - f!o_sent; 
nodeCnd] • d©ma_supl tt 1 m] + f lo_sent ; 
nodeCnd] . requ^avai Cd_t im3 + flo_sent; 


node [ nd 3 . dema__supl C 1 1 m3 : = 
nodeCnd] . requ_avai [d_t im] 2 = 
WRITELN (’tim d ’ , 


d_t tau 2 * 

END; 

WITH spnC } store 3 DO 
BEGIN 

sup := sup - f 1 o_sent 
FOR kk := 1 TO dnn DC 
ass ignt j store , kk3 : ~ 


ass ignC jstore,kk3 - flo_sent; 


WRITELN (’tim s * ,s_tim:2>i 
nodeCnd 3 . dema_supl 1 1 im3 : ~ 


nodeCnd] - dema_supl tt im3 "* f lo_ sen ^ ; 



THEN nodeCnd] . requ avaiCs timj := 

odeCnd 3 .roqu_avaiUs_.ttm 3 - fio sent; ~ ~ 

END; 

(WRITE arcs aray ok fwdl 
kk ;= d_tim,- 
rk := kk; 

FOR j j : = no DOWNTO 1 DO 
WITH arcC arcs_aray [ j ) 3 3DO 
BEG-IN 

FOR tt := kk DOWNTO 

kk-round(di st> DO 


res_capC 1 1 3 -f 1 o_sent ; 
loadCtt} + flo sent; 


arcs_arayC ) } 3 : 2 ) ; 


res_capt 1 1 3 : 5 « 0 , ' Load 


BEGIN 

res_cap[tt] : = 


loadCtt3 : = 

I 


1 , 1 oad[ 1 1 3 » J : 0 ) j 


END; 


END; 

WRITELN (’Arc no 

j 

I 

WRITELN (’Res capa # t 

WRITE ( 'rk : * , rk: 5:G> ; j 

rk ;= rk - dist; I 

WR I TELN ( rk : 3 : 0 > ; 
kk : = round! rk) ; 


WRITELN! 'flo sent ’ , f lo_sent : 4 : 0 > ; 
END C ass > 0.001 > 

ELSE GOTO 250;! ass <= 0.001 } 


END 

ELSE GOTO 250; (node 2 not in dem_nod_set) 

node!s_nod] . cum_dem_sup_g := 

nodeCs_nod3 . cum_dem_sup_g - f lo_sent ; 

FOR J J := 1 TO no DO 
WITH arc[arcs_aray[ j) 33 DO 
BEGIN 

tmp_capa := tmp_capa - flo_sent; 
IF tmp_capa < 0.001 
THEN deleted := TRUE; 
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node [ node 2 3 . cum_dem_sup_g := 

node C node 2 J - cum_dem_sup_g + flo_sent; 

IF node[node2] . cumdem supq > -0.001 
IHEW dem_nod_set : = dem_nod_set - Cnode2 ] ; 
WRITELN ( 'flag 2 • ) j 
GOTO 100; 

END; CIF ft in arc_flo_st & not arcC f t ] . de le td> 
ENO; Cfor ii ;= fpoint.. & WITH arct i l 3 do} 

WRITELN; 

Cs_nod had surplus suplies which could not be absorbed 
WITH current assignments} 

IF s_nod <= no_nod 
THEN IF no > 0 
THEN BEGIN 

arc[arcs_aray [no] 3 . enter := FALSE; 
min_arc_capa ; = prev_mi n_arc_cap ; 
cur_nod := arc[arcs_aray [no] 3 . nodel ; 
no : = no - 1 ; 

WRITELN ('Backtracking'); 

GOTO 75; 

END 

ELSE BEGIN 

s nod : = s_nod + 1 ; 

WRITELN ('Futile search, 

s_nod Quit,, surplus assumed') 

END 

ELSE BEGIN C s_nod > no_nod > 

WRITELN ('Failure with current assignments, 

discarding them* > i 

s_nod s - 1 ; 

FOR ii : = 1 TO snn DO 

FOR jj := 1 T0 dnn 00 

IF (spnCn3.sup50.001) AND (dmnC JJ 3 . dem>0 . 001 > 

THEN assignC i i . H 3 ;= smaller 

(spnC i i ] . sup.dmnC j j 3 .dem) ; 

END; 

cn ;* cn + li 
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100 ; 


350 : 


375 : 


cn t = cn + 1 ; 

UNTIL (flo_sent > 0.001) OR (cn>400>: 

WRITELN ( * f lag5 # ) : 
rk : = s_t i m ? 

IF f 1 o_sent > 0.001 
THEN BEGIN 

WRITE ( out 2 , f 1 o_sent : 6 : 2 , * *); 

FOR i l : = 1 TO no DO 

WITH arc[arcs__aray[ i i ]] DO 

BEGIN 

WRI TE (out2 , nodel s 8 , rk : 8 : 2) : 
rk : = rk + cf i s t ; 

END? 

IF no <> 0 

THEN WITH arc[arcs_aray[no}] DO 

WRITELN (out2,node2:8,rk:8:2>: 

END; 

END C while dem__nod_s e t <> [ ] } 

ELSE C mincoat infeasible > 

WHILE (dem_nod_set <> []) AND (cn<180) DO 
BEGIN 

WRITE ( 'Dem nod set : ’); 

FOR it : = 1 TO no_nod DO 

IF ii IN dem_nod__set 
THEN WRITE <ii:2>; WRITELN; 

FOR l i ;= 1 TO snn DO 
FOR j ) i= 1 TO dnn DO 

IF (ass l gn[ i i f j j ] > 0) and (dmn[jj].nd IN dem__nod_set ) 
THEN GOTO 350; 
s_nod := spnCiil.nd; 
d__nod : = dmntjjl.nd; 

avai l_t lm := dmntjjl.tim - spnCiil.tim; 

f 1 o_sent := smaller (spn[ l l ] .sup, dmnC j J 3 . dem) ; 

flo_sent := smaller (f lo^sent, ass i gnt i l , j j ] > ; 

istore := it? jstore s= jj; 

rk : = dmn[ j j 3 . t im; 

excl^arc i = C 3 « tried := FALSE: 

f i nd_short_ex__a ray ( s_nod , d__nod ,exci__arc ,arcs_aray , no) t 
C Here arcs^aray's storage is in reverse order > 
mot ion t i m : = 0 ; 
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380 : 


FOR n t= 1 TO no DO 

motion^t im : = motion_tim + ar c [ arc s_aray [ i i 3 3 . d i s t * 

I F (mot ion_tim > aval l_t im) AND NOT tried 
THEN BEGIN 

WRITELN (’Time infeasible’); 
tried := TRUE; 

IF excl^arc <> [3 
THEN BEGIN 

e x c 1 __a re : = C 3 t 
goto 373 ; 

END; 

END; 

WR I TELN < * s_nod * , s_nod : 3 , ’ d_nod # , d_nod : 3 ) ; 

WRITE (’Short aray ’ ) ; 

FOR it := no DOWNTO 1 DO 

WRITE (arcs__aray[ i i 3 : 4) ; 

WRITELN; 

C We check the residual capacities of the arcs, and time f eas i 3 
travl_tim := 0; 

FOR i 1 := 1 TO no DO 

WITH arc[arcs_aray[ i i 3 3 DO 

BEGIN 

t i m_d t := rk~ t rav 1 _t i m; 
t i m_s t : = t i m_d t - diet; 

IF tim_st < 0 THEN tim_st r = 0; 
passed ; = TRUE; 

FOR tt t ** round ( 1 1 m_d t ) DOWNTO round ( t i m_B t > DO 
BEGIN 

IF res_cap[tt3 <“ 0 
THEN BEGIN 

t im_dt := tt - 1 ; 
tim_st : = tirn_dt - dist: 
passed := FALSE; 

END 

ELSE f 1 o_sent : = sma 1 i e r ( f 1 o_sent , res_capC t t 3 > f 

END; 

IF NOT passed 

THEN IF degree [arc Car cs__arayC i i 3 3 . node 2 3 > 1 

THEN e x c I _a r c := excl__arc + [arcs^aray [ i i 3 3 
ELSE tried : = TRUE; 
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tra victim : = rk - 1 1 m_s t j 

IF NOT passed AND ( ( travl_t im<avai I_t im) OR tried) 
THEN GOTO >80 i 

s i ack_t i m := aval i_t lm - travl_t imi 
IF (slack__t im < 0 > AND (mot ion__t im<avai l_tim) 

AND NOT tried 

THEN BEGIN 

WRITE < ' Exc 1 Arc ’ ) ; 

FOR j j : = 1 TO 2*no_arc DO 
I F i j IN exc 1 _a r c 
THEN WRITE ( j j : 4> ; 

WRITELN; 

WRITELN (’Time infeasible (temp)’); 

GOTO 373; 

END; 

C else if mo t i on_i i m > tim then delays ) 

END; 

IF aval l_t im < travl__tim 
THEN BEGIN 

WRITE (’Time i n f eas i b i 1 i t y ’>; 
delay ; = t rav 1 _t l m - a v a s 1 __ t i m ; 

END 

ELSE delay : = 0; 

WRITELN (’delay ’ t de 1 ay : ? : 2 > ; 

Updat e__de 1 ay_i nf o ( de 1 ay , f i o_sent , de 1 ay_f req, wt d__de lay ,ptr) ; 
FOR i i * : = no DOWNTO 1 DO 
WITH arc[arcs_aray[ i i ]] DO 
BEGIN 

FOR tt ; = r ound ( t i m_s t ) TO round ( t i m__dt > DO 
BEGIN 

res_cap[ 1 1 ] : = re s__cap[ 1 1 ] - f 1 o_sent ; 

loadftt] ! = loadCtt] F flo^sent; 

END ; 

WRITELN (’arc no ’ , arcs_aray C i i ] r 2 , 

’ Res capa * f res^capC 1 1 T : 3 : 0 , load * , 1 oadC 1 1 3 : 3 ; 0 ) ; 

WRITE ( * t i m s t i • * tim_sts 3s0> « 

WRITELN (t im_dt:>:0) ; 

’ per i od_cos t : = per i od_cost + f 1 o_sent* cos t_g ; 


END 
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WRITELN <*fio sent * , f i o_sent : 4 r 1 > ; 

assignfistore, )store]:= assign[ istore, jstore] -f 1 o_sent t 
spnC i store 3 . sup : = spn[ i s t ore 3 . sup - f 1 o_ sent; 
time := spn[ l s tore ] . t ,im; 

WITH node [ s_/iod 3 DO 
BEGIN 

dema__sup 1 [ t i me ] := dema_sup 1 C t i me 3 - f 1 o__sent ; 

r e qu ava i [ t i me 3 := requ__ava i [ t l me ] - f 1 absent ; 

cum_dem_sup_g : = cum_dem__sup_g - f 1 o_sent ; 

END: 

dmn C j s tore 3 . dem := dmn[ j store 3- dem - fio__sent; 
time := dmnC j s tore 3 . t im: 

WITH node[d_nod3 DO 
BEGIN 

dema^sup 1 [ 1 1 me 3 s = dema_9up 1 [ 1 1 me 3 + f 1 o_sent : 
requ_ava i [ t i me 3 := requ^ava i [ 1 1 me 3 F f 1 o_sent ; 
cum „^ em __ s U P__9 : - cum_dem_sup_g + f I o_sent ; 

IF cum_dem__9up_g > -0.001 

THEN dem_nod_s e t .*= dem_nod_s e t - [d_nod3; 

END? 

IF f 1 absent > 0.001 
THEN BEGIN 

WRITE ( ou 1 2 , f 1 o__sent : 6 r 2 , * *>: 

FOR i i := no DOWNTO 1 DO 
WITH arc[arcs_aray[ i l 33 DO 

WRITE (out2,nodel :8, tim_st :8:2>i 
IF no <> 0 

THEN WITH arc tarcs_aray [ 1 3 3 DO 

WR I TELN( out 2 , node 2 : 8. < t im__s t+d i s t > : 8 ; 2 > 

ELSE WRITELN ( out2 , s_nod : 8 , spnC i s tore ] . t i m: 6 . d_nod * 8 , 
dmnC jstore 3. tim:6>: 


END; 

cn j = cn +1 ; 


END; C while <dem_nod_set <> C3> and (cn<180> do! 

IF dem_nod_set <> C3 THEN 

WRITELN (’Exiting while loop !! goods demand unsat i sf i ed ’> z 
WRITE (outJ,period_no:4 f periodic os t : 12 : 2 ) ; 
total cost := total_cost + period_cost; 


500 : 
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C Empty wagon scheduling 

y 

WRITELN; WRITELN < * Empty wagon schedul ing’ > ; 
WRITELN (’Period no : ' ,period__no.*2) ; 

per l od_cost := 0; 
d em __ n °d__ s e t := C ] ; 

FOR li := 1 TO no__nod DO 

WITH nodeC i i 3 00 

BEGIN 

IF penod^no = 1 
THEN BEGIN 


cum_dem_sup__w := 0; 

FOR tt := 1 TO peri od__i en th DO 
IF requ_avai [tt] > 0.001 

THEN tum_dem_sup_w := cum_dem_sup_w + requ_ava i C t t 3 

END 

ELSE C period^no > 1 } 

FOR tt := (per iod_no-l )*period_ienth + 1 
to period_no*period_lenth DO 
i Consider demand for wagons in current time period > 

BEGIN 

IF requ__ava i [ 1 1 3 < -0.001 
THEN BEGIN 


+ requ^avai C tt 3 ; 


cum__dem_sup__w : = cum_dem_sup_w 

dem_nod_set : = dem_nod_set + Cii3; 
END; 


t Consider suply in wagons in previous time period > 
IF requ__avai [tt ] > 0.001 
THEN cum_dem_sup_w : = c u m_d e m_s u p__w 


+ requ_avai Ctt 3 ; 


END; 


END; 

IF period^no = 1 
THEN BEGIN 

WRITELN (’Empty wagon scheduling is assumed 
to be done in the previous time horizon'); 

WRITELN (out?); 


{ Also mention the wagons for which assume this > 
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GOTO 600; 

END; 

WRITELN < f Node cum_dem_sup__w' ) ; 

FOR ii : = 1 TO no__nod DO 

WRITELN <ii:4 t nodeC nl. cum_dem__sup_w : 8:1); 

WRITE < * Dem nod set ; 1 > ; 

FOR ii := 1 TO no^nod DO 
IF ii IN dem_nod_s e t 
THEN WRITE <ii:2>; WRITELN; 
dnn := 0; i No. of dem nodes > 
snn := 0; C No. of sup nodes > 

FOR ii : = I TO no_nod+I DO 
BEGIN 

ddegC i l ] : = 0 ; 

sdegC l i 3 : = 0 ; 

END; 

FOR ii := 1 TO no_nod DO 
WITH nodeC i i3 DO 
BEGIN 

FOR tt := per i od_i enth + 1 TO T__ DO 
IF requ_avai[tt] < -0.001 
THEN BEGIN 

dnn : = dnn + 1 ; 

ddegC i l ] := ddegC l i ] + 1; 

WITH dmnCdnn] DO 
BEGIN nd : = ii; 

dem := abs< requ__ava l C tt 3 ) ; 
t im : = t t ; 

END; 

END 

END; 

FOR ii := 1 TO no_nod DO 

WITH nodeC i i 3 DO 

BEGIN 

FOR tt : = 1 TO T_ DO 
IF requ^avai C tt 3 > 0.001 
THEN BEGIN 


snn : = snn + 1 ; 

sdegC i l 3 sdegC i i 3 + 1; 



550: 


575 : 
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WRITE ( dmnC 1 1 3 . dem: 5:0); 

WRITELN; 

WRITE ( 9 9 > ; 

FOR li := 1 TO dnn DO 

WRITE <dmnCi i3.tim:3>; 

WRITELN; 

FOR s_nod : = 1 TO no_nod DO 

FOR )j:= spt[s___nod] TO spt [ s__nodr 1 ] - 1 DO 

BEGIN 

WR I TE < spnC j } ] . nd : 5 , spn[ j j ] . sup : 3 : 0 t spnC J J 3 . 1 1 m : 5 ) ; 

FOR cur__nod := 1 TO no__nod DO 

FOR kk:= dpt[cur_nod] TO dpt [ cur__nod+l ] - 1 DO 
WRITE (assign[jj,kk3:3:0>; 

CWR I TE ( short est_d i st ( spnC j j 3 . nd , dmn[kk 3 . nd> : 3 : Q > ; > 

WRITELN; 

END; 

WRITELN (out2 f 'Empty Wagon Movement Period no : # , per i od_no : 2 ) 

WRITELN (out2 t *Flow Node Time *>; 

CH : = 1 ; s^nod : = 1 ; 

WHILE (dem_nod_set <> C3> AND (cn<200) DO 
BEGIN 

WRITE (’Dem nod set : ’>; 

FOR li ; = 1 TO no_nod DO 

IF li IN dem_nod_se t 
THEN WRITE <ii:2>; WRITELN; 

FOR ii : = 1 TO snn DO 
FOR j j := 1 TO dnn DO 

IF ass i gnC i i , j j 3 > 0.001 
THEN GOTO 550; 
s_nod := spnCii3.nd; 
d_nod := dmnCj}3.nd; 

aval l_t lm : = dmnCjjl.tim - spnCii3*tim; 

flo_sent ; = smaller < spnC i i 3 . sup , dmnC j j 3 . dem) ; 

flo_sent := smaller (f lo^sent , ass l gnC i i , j j 3 > ; 

istore := n; jstore t = jj; 

rk : = spnC i l 3 . t im; 

excl__arc : = [3; tried := FALSE; 

find short_ex_aray < s_nod t d_nod » exc l_arc , arcs^aray , no) ; 


C Here arcs_aray’s storage is in reverse order > 
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mo 1 i on__t i m : = 0; 

FOR 11 := 1 TO no DO 

mo 1 1 on_t i m := motion^tim * arc Carcs^a ray C l i ] ] . d i s t ; 
IF <mot ion_t lm > aval AND NOT tried 
THEN BEGIN 

WRITELN < ’Time infeasible*); 
tried := TRUE; 

IF exc I_arc <> [] 

THEN BEGIN 


exc i_ar c : = C 3 ; 
goto 575; 

END; 

END; 

WRITELN < ' s_nod * , s_nod : 5 , * d__nod * , d_nod : 5 > ; 
WRITE ( 'Short aray * > ; 

FOR n : = no DOWNTO I DO 

WRITE (arcs__aray[ l l ] : 4) ; 

WRITELN ; 

t We check the residual capacities of the arcs, 
and time feasibility before sending the flo } 

tra victim : = 0; 

FOR li := no DOWNTO 1 DO 
WITH arc[arcs_aray C i l ] ] DO 
BEGIN 


580 : 


t i m_s t := rk+travi_t im; 
passed : = TRUE; 

FOR tt := round ( 1 1 m__s t > TO round ( 1 1 m__s t+d l st > DO 
BEGIN 

IF res__cap[tt] <= 0 
THEN BEGIN 

t l m__s t : - 1 1 + 1 ; 
passed : = FALSE; 

END 

ELSE f 1 o_sent := sma 1 1 er < f 1 o_sen t , res_capC 1 1 ] ) 

END; 

IF NOT passed 

THEN exc l_arc := excl_arc + Car cs_aray C i i 3 3 ; 
tra victim : = (tim_st-rk> + dist; 

IF NOT passed AND ( ( t rav l_t im<ava l l_t lm) OR tried) 


THEN GOTO 580; 
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AND NOT tried 


END; 


slack_t lm := aval l_t im - t rav I__t i m ; 

IF (slack__tim < 0> AND (mot ion_t inKavai 

THEN BEGIN 

WRI TELN ( ’Time 1 nfeas 1 bi e (temp)*); 
GOTO 575; 

END; 

C else if mot 1 on_t 1 m > 1 1 m then de 1 ays } 


IF aval l^tim < travl__tim 
THEN BEGIN 

WRITE ('Time l nfeas l b i 1 l ty * ) ; 
delay := trav 1 __t i m - aval l_t im; 

END 

ELSE .de lay ; = 0 ; 

WRITELN (’delay * , de 1 ay : 5 : 2 ) ; 

Updat e__de lay_i n f o ( de 1 ay , f 1 o__sen t , de 1 ay__f req , wt d__de lay , pt r ) 
FOR i l := no DOWNTO 1 DO 
WITH arc[arcs_arayt i i 13 DO 
BEGIN 

FOR tt := round ( 1 1 m__s t ) TO round ( 1 1 m_s t+d i st ) DO 
BEGIN 

res_cap[ 1 1 ] : = r e s__cap [ 1 1 ] - f 1 o_sent ; 

loadCtt] : = loadttt] + f lo__sent j 

END; 

WRITELN ('arc no * , arcs^aray C i i 3 : 2 , 

* Res capa ' ,res__cap[tt]:):0, * load ' # 1 oadC t t ] : 3 : 0 > ; 

WRITE ( ' t im s t : ' , t im_s t : 3 : 0 ) ; 

t 1 t : = t i m_s t + d l s t ; 

WRITELN <tim_dt:3:0) ; 

per iod_cost : = per iod_cost + f 1 o__sen t*cos t_w ; 

END; 

WRITELN (’flo sent ’ , f 1 o_sent : 3 : 0 ) ; 

ass i gnC i s tore , j s t ore 3 : - ass i gnC i store , j store 3 ~flo_sent; 
spn[ istore] .sup := spn[ i s tore 3 . sup - flo_sent; 
time : = spnC istore] . t im; 

WITH nodeCs^nod] DO 
BEGIN 


requ^ava i [ t l me ] := r equ_avai [ t ime 3 - flo_sent; 
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cum_dem_sup_w : = cum_dem_sup_w - f 1 o_sent ; 

END; 

dmnC j store 3 . dem := dmn[ ) store ]. dem - fio_sent; 
time := dmnC j s tore 3 . t im'; 

WITH node£d_nod] DO 
BEGIN 

requ__avai [t ime] : = requ_avai Ct ime] + f 1 o_sent ; 
cum__dem_sup__w : = cum_dem__sup__w + flo_sent; 

IF cum__dem_sup_w > - 0.001 

THEN dem__nod_set : = dem_nod_set - C d___nod J ; 

END; 

IF f 1 o^sent > 0.001 
THEN BEGIN 

WRITE (out2,f lo_sent:6:2, * ’>; 

FOR li := no DOWNTO 1 DO 
WITH arcCarcs_aray [ l i 3 3 DO 

WR ITE <out2,nodel:8,ti m__s t : 8 : 2 ) ; 

IF no <> 0 

THEN WITH arc Care s_aray C 111 DO 

WRITELN<out2,node2:8, < t im_st+di st > :8:2) 
ELSE WRITELN ( out 2 , s_nod : 8 , spnC l store 3 . t im: 6 f 
d_nod : 8 , dmnC jstore3.tim:6> ; 

END; 

cn : = cn +1 ; 

END; £ while (dem_nod_set <> C3> and (cn<200) do> 

IF dem_nod_set <> [ 3 

THEN WRITELN (’Exiting while loop !! emp wag dem unsatisfied* ) ; 
WRITELN ( out 5 , per i od_cost : 1 2 : 2 > ; 
totai__cost : = total__cost + period_cost; 

600: IF period_no* period_ienth < T_hor 

THEN BEGIN 

period_no : = period_no + 1; 

GOTO ?0; 

ENO; 

REWRITE (outl , 'outdelay.dat ’> ; 

WRITELN < outl , 'Delay Frequency % Delay wtd % Delay'); 
cumsum := 0; tot_delay := 0; 

FOR ii ;= 0 TO ptr DO 


BEGIN 
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cumsum : = cumsum + de iay_f reqC i 1 1 ; 
tot_de 1 ay : = to t__de lay + wtd_de i ay C l i 3 ; 

END; 

FOR ii := 0 TO ptr DO 

WRITELN < out 1 , l i : 4 • ’ * , de iay_f reqC l i ] : 4 f delay 

cumsum*! 00 : 1 2 : 2 f 


wtd_de lay C i i 1 / tot_de 1 ay*l 00:12:2); 
WRITELN ( out 1 ) ; 

WRITELN (outl , ’Total Frequency 
WRITELN ( outl t ’Total Wtd delay 
WRITELN (outl , ’Average Delay 
CLOSE (outl); 

CLOSE ( ou 1 2 ) ; 

REWRITE (out) , ’outload.dat ’ ) ; 
WRITELN <out3L’Arc no Load 
FOR ii : = 1 TO 2*no___ar c DO 
WITH arct l l ] DO 
BEGIN 


* , cumsum ; 4 ) ; 

’ , t ot_de 1 ay : 6 : 2 ) ; 

* , tot_de lay/cumsum 


Time'); 


max load := 0; 

FOR tt ;= 1 TO T_hor DO 
BEGIN 

IF maxload < loadCtt] 

THEN BEGIN 

max load : = loadCtt] ; 
t ime : = t t ; 

END; 

END; 

WRITELN ( out 5 , l l : 6 .max load :8:0*time:8); 

END; 

CLOSE ( out 3 ) ; 

REWRITE(out4, 1 da l 1 y 1 oad . d ’ ) ; 
per l od_no : = 1 ; 

WHILE peri od_no * per 10 d__ 1 e n t h <= T_hor DO 
BEGIN 

WRITELN (out4 . ' Per l od__no : f , per i od__no : 2 ) 
WRITELN ( out 4 , ’ Arc No. Time’); 

WRITE (out4, * f ) ; 

FOR tt : = (per iod_no-l )*per iod_lenth + 1 
peri od__no * period_lenth 


f reqC i i 1 / 


: 6 : 2 ) ; 


to 

DO 


WRITE < out 4 , 1 1 s 4 ) ; 
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WRITELN <out4>; 

FOR 11 := 1 TO 2 * no_arc DO 

WITH arcC l i ] DO 

BEGIN 

WRITE ( out 4 ,ii:4, ' ’); 

FOR tt := ( per l od_no-l > *per l od_l enth +1 to 

per iod_no * period_lenth DO 
WRITE <out4, loadCtt] :4: 1 ) ; 

WRITELN <out4> ; 

END; 

penod_no := per iod_no + l- 

END; 

CLOSE (out4> ; 

WRITELN (out?); 

WRITELN (out5, 'Total Cost : ' , tota I_cost : 1 2 : 2 > ; 

CLOSE (out?); 

END. 



APPENDIX C 


PROBLEMS SOLVED AND THEIR RESULTS 

Here we present results of the 6 more problems. All the 
problems solved, took a CPU time in the range of 27-32 seconds of 
HP 9000. The results obtained are more less similar to those 
presented in the chapter 3. However here we have noticed that for 
the problems having tighter capacities of the arcs have not only 
more delays but also more costs of transportation, which appears 
clearly in graphs. 


PROBLEMS 3 AND 4 - 

The problems 1 & 2 are solved as 3 (instead of 4) period problems. 
Time horizon = 96 time units. Period Length = 20 time units 
FINAL RESULTS FOR THE PROBLEMS 3 AND 4 - 

Table C. 1 Final Results for the Problem 3 and 4 


Problem 3 


Empty klagons 

Transportation Costs 

Delays 

128 

36,446 

4.26 

153 

35,968 

3.97 

178 

35,565 

3.85 

203 

35. 105 

3.66 

128 

39,677 

4.43 

153 

38,570 

4.12 

178 

38,261 

3.94 

203 

37,547 

3.73 


Problem 4 






0 


INPUT FOR THE PROBLEMS 5 AND 6 
Table C. 2. 1 Supply at various nodes 


Supply Node 

Supply, time (day) 

5 Kanpur 

4,5; 6,4; 7,5; 5,6; 10,7; 5,9; 4,10; 

5 Calcutta 

6.1; 7,4; 12,5; 5,6; 2,11; 8,12; 

7 V i shakhap . 

8.2; 5,5; 4,2; 6.5; 7,6; 5,9; 6.10; 

8 Madras 

5,1 ; 4,2; 8,7; 5,8; 4,11 ; 6.12; 

1 5 Ahmadabad 

2,1 ; 8.5; 5,7; 6,8; 6,11; 7,12; 

1 4 Manmad 

5,1; 4,2; 5,5; 4.6; 5.9; 4,10; 

16 Nagpur 

8.4; 5,1; 5,7; 4,8; 7,9; 5,11; 10,12 

19 Bhopal 

7,1; 5.2; 10,5; 4,5; 6,6; 8,9; 5,10; 

20 Indore 

5,2; 6,1; 2,7; 8,8; 8,9; 5,10; 4,11; j 


Table C. 2. 2 Demands at various nodes 


Demand node 

Demand, time (day) 

2 Jhansi 

4,5; 6,4; 4,5; 6,6; 8.9; 5,10; 

18 Hydrabad 

5.2; 6,1; 5.7; 6.8; 6,11; 7,12; 5,9; 6,10 

4 Lucknow 

8,1; 7,4; 8.5; 6,5; 7,6; 8,7; 5,8; 4,11; 


6,12; 5,9; 4,10; 

9 Bangalore 

8,4; 5,1; 4,2; 5,5; 4,6; 7,9; 5,11; 10,12 

6 Jamshedpur 

5.1; 4,2; 5,7; 4.8; 2,11; 8,12 

1 Delhi 

5,1; 4,2; 8.5; 5.6; 4,5; 2,7; 8,8; 8,9 


5,10; 4,11 


1 2 Bombay 


8,2? 15.5; 7,5; 5,6; 10,7; 5.9; 4,10 














FINAL RESULTS FOR PROBLEMS 5 AND 6 

Table C. 2. 4 Final Results for the Problem 5 


Empty Wagons 

Transportation Costs 

Delays 

205 

22,471 

4.56 

250 

21 . 568 

CD 

04 

• 

255 

21 ,042 

4.15 

280 

21,042 

4.17 


Table C. 2. 5 Final Results for the Problem 6 


Empty Wagons 

Transportation Costs 

Delays 

205 

22,595 

5.64 

250 

21,480 

5.64 

255 

21,155 

5.55 

280 

21,155 

5.52 
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Transport. Cost Vs. Wagons Available 



Xol of Wnjoni Att liable 

FIG. C.3 

Delays Vs. No. of Wagons Available 
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Ko. of Wif ora A-rtUiils 


FIG. C.4- 
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INPUT FOR THE PROBLEMS 5 AND 6 
Table C. 2. 1 Supply at var l ous nodes 


Supply Node 

Supply, time (day) 

5 Kanpur 

4.3; 6,4; 7,5; 5,6; 10,7; 5,9; 4,10; 

5 Calcutta 

6,1; 7,4; 12,5; 5,6; 2,11; 8,12; 

7 Vishakhap. 

8,2; 5,3; 4,2; 6.5; 7,6; 5.9; 6.10; 

8 Madras 

5,1; 4,2; 8,7; 5.8; 4,11; 6,12; 

1 5 Ahmadabad 

2,1; 8,3; 5,7; 6,8; 6,11 ; 7,12; 

1 4 Manmad 

3,1; 4,2; 3.5; 4.6; 3.9; 4,10; 

16 Nagpur 

8,4; 5.1; 5,7; 4.8; 7.9; 5.11; 10,12 

19 Bhopal 

7,1; 5.2; 10,3; 4.5; 6,6; 8,9; 5,10; 

20 Indore 

5,2; 6,1; 2,7; 8,8; 8,9; 5,10; 4,11; 


Table C. 2. 2 Demands at various nodes 


Oemand node 

Demand, time (day) 

2 Jhans i 

4,3; 6,4; 4,5; 6.6; 8,9; 5,10; 

18 Hydrabad 

5,2; 6,1; 5.7; 6.8; 6,11; 7,12; 5,9; 6,10 

4 Lucknow 

8,1; 7,4; 8,3; 6,5; 7,6; 8,7; 5.8; 4,11; 


6,12; 3,9; 4,10; 

9 Bangalore 

8,4; 5,1; 4,2; 3,5; 4,6; 7.9; 5,11; 10,12 

6 Jamshedpur 

3,1; 4.2; 5,7; 4,8; 2,11; 8,12 

1 Delhi 

5,1; 4,2; 8,5; 5.6; 4,5; 2.7; 8,8; 8.9 


1 2 Bombay 


5.10; 4,11 

8,2; 15.3; 7.5; 5.6; 10.7; 5,9; 4.10 












FINAL RESULTS FOR PROBLEMS 5 AND 6 

Table C. 2. 4t Final Results for the Problem 5 



Table C.2«5 Final Results for the Problem 6 


Empty Nagons 

Transportation Costs 

Delays 

205 

22,393 

5,64 

250 

21,480 

5.64 

255 

21,153 

5.33 

280 

21,153 

5.32 
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INPUT FOR THE PROBLEMS 7 AND 8 

Table C. 3.1 Supplies at various nodes 


Supply Node 

Supply, time (day) 

1? Bhusaval 

6,1; 7.2; 1,6; 4.7; 6,11; 7,12 

20 Indore 

8,?; 5,4; 8,8; 7,9; 5,10; 8,12; 5,15 

? Kanpur 

4,5; 5,1; 4.6; 5,7; 4,8; 1,14; 4,15 

4 Lucknow 

4,2; 2.5; 4,9; 2,10; 8,11; 7,12; 5,15 

9 Bangalore 

4,4; 5,5; 4.1; 4,6; 5,7; 4,14; 3,16; 4,11 

6 Jamshedpur 

8,2; 7,3; 5,4; 8,8; 5,9; 4,12; 2,1? 

1 Delhi 

1,5; 4,1; 6,10; 7.6; 4,14; 5.15 


Table C. 3. 2 Demands at various node 


Demand Node 

Demand, time (day) 

19 Bhopal 

6,2; 4,10; 4,6; 5,15; 7,11 

1 3 Ahmdabad 

1,5; 4,1; 3,8; 5,9; 4,13; 4,14 

7 Vishakha. 

6,4; 1,10; 2,6; 5,7; 3,11; 5,12 

1 1 Pune 

1.2; 2,3; 3.7; 4.8; 1,9; 2,14; 6,15 

17 Balharsha 

7,5; 2,1; 5.10; 3,6; 8,13 

1 0 Daund 

1.4; 8,9; 5.11; 3,12 

1 4 Manmad 

3.2; 3,3; 2,7; 6.8; 3,13; 1.15 

16 Nagpur 

4,5; 5,1; 5,10; 5,6; 1,15; 2,11; 5,12 

2 Jhans i 

1,3; 7,4; 4,8; 4,9; 3,13; 9,14 

18 Hydrabad 

3,1; 4,2; 5,6; 7,7; 4,11; 4,12 


The time horizon = 1? days and period length = 5 days. 
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Table C. 3. 3 Table for capacities of the arcs. 


Arc 

No . 

n 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

Pr. 

7 

B 

5 

25 

8 

20 

24 

5 

21 

6 

18 

18 

7 

6 

23 

Pr. 

8 


5 

20 

8 

20 

20 

5 

20 

6 

15 

10 

6 

7 

22 



FINAL RESULTS FOR THE PROBLEMS 7 AND 8 
Table C. 3. 4 Final Results for problem 7 


Empty Wagons 

Transportation Costs 

De lay 

157 

21,753 

5.25 

182 

20.790 

5.03 

207 

19,786 

4.61 

232 

17,843 

4.62 


Table C.3. 4 Final Results for problem 8 


Empty Wagons 

Transportation Costs 

Delay 

157 

22,187 

6.06 

182 

21,241 

5.63 

207 

20,223 

5.32 

232 

18,300 

4.72 























DISCUSSION 


Here we give minutes of discussion during the thesis presentat 1 on . 
Q. What is the smallest unit of time for which scheduling is 
done ? What is the basis of deciding the its level ? 

Ans. We have used 1 unit time = 4 hrs. in our scheduling problem. 
The criteria for deciding the smallest unit of time are 
1. Travelling time for arc with smallest dist Z. Smaller the unit 
time finer the level of scheduling and more computational time. 

Q. Where do we use simulation in the methodology ? 

Ans. During routing we try incorporate to have essence of real 
life logic and relationships. Further, we perform experiments we 
the model developed so as to decide what fraction of capacity of 
arcs per period should be taken as the capacity of arc per unit 
time. 

Q. What is alternative methodology that can be used ? 

Ans. One can try to use Lagrangian Relaxation by relaxing some 
constraints, dualizing them and exploit some standard structure 
out of the problem. Also some other ways of dividing the problem 
into smaller subproblems, methods like Benders* Decomposition may 
be tried. Also one can have analogy between our problem and job 
shop scheduling problem. Wagons are analogous to the jobs. 
Wagons moving over arcs is similar to the jobs loaded on the 



machines. Travelling time being analogous to processing time of 
jobs. Delivery dates for the goods are similar to due dates of 
the jobs. However it should be noted that scheduling many jobs on 
more machines <more than 5 > is itself a tough problem to solve. 

Q* How the costs are expressed? Can they be expressed in terms 
of cost of moving a wagon per arc? 

Ans. Costs are expressed as cost of moving a train of wagons per 
arc. It can be expressed in terms of wagons also for going into 
finer details of scheduling. It will involve some component of 
cost which is fixed and some component which is proportional to 
the number of wagons in train. 

Q. Is it necessary to take cost proportional to the distance ? 
Ans. It is not necessary to take cost proportional to the 
distance. In real world also cost may vary not only with distance 
but also depending on whether the route is electrified or uses 
diesel or coal engines. The algorithm will work such input also 
without any further modification. 



